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Abstract 

We analyze two models of random geometries : planar hyper-cubic random surfaces 
and four dimensional simplicial quantum gravity. 

We show for the hyper-cubic random surface model that a geometrical con- 
straint does not change the critical properties of the model compared to the model 
without this constraint. We analyze the phase diagram for the model with extrinsic 
curvature. 

For four dimensional simplicial quantum gravity we find that in the large volume 
limit the leading contribution to the entropy does not depend on the underlying 
topology. We find for the first time a strong back-reaction of the geometric sector 
for matter coupled to gravity. 
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Chapter 1 

Introduction 



At the beginning of this century two theories forced a new understanding of physics. 
Quantum theory on the one hand describes the physics on small scales, where 
"small" scales do not only appear on the atomic level but also for macroscopic 
objects, such as white dwarfs or super-fluid Helium. On the other hand we have 
general relativity describing large distance physics which is dominated by the grav- 
itational field. 

As of today all attempts to unify these formalisms failed. This is very unfor- 
tunate because with such a theory at hand it may be possible to overcome a weak 
point of general relativity: it predicts space-time singularities and at the same time 
breaks down at these singularities. In particular, the Big Bang could have produced 
anything, as far as classical general relativity is concerned. A theory of quantum 
gravity on the other hand might be able to produce some boundary conditions that 
can tell us why the universe looks like it does. 

It is however very difficult to construct such a theory. The scale in quan- 
tum gravity, where quantum effects become important, is set by the Planck mass 
Mp ~ 10^^ GeV. Typical energies in particle physics experiments are, however, sev- 
eral orders of magnitude smaller, which is why no experimental results for quantum 
effects of gravity are available. Therefore one has no no guidance from experiments, 
which, for example, in the case of the standard model provided very important in- 
sight necessary to construct the theory. 

String theories are one attempt to unify the fundamental interactions in nature 
(see for example These theories avoid the singularities which usually plague 
quantum field theories by considering the propagation of strings rather than point 
like particles. Differently from point particles which propagate along paths strings 
propagate along the world-sheet, a two dimensional surface. Is is very interesting 
that string theories contain spin 2 particles, i.e. gravitons, which makes them a 
promising candidate for quantum gravity. 

General Relativity and string theory have in common that they both are in a 
certain sense geometric theories. General relativity describes gravity as a curved 
four dimensional space time metric, strings propagate along two dimensional sur- 
faces. As we will describe in more detail below quantum gravity and quantum 
string theory both can be looked at as a statistical system of geometries, random 
geometries. Such systems are, unfortunately, very difficult to deal with analytically, 
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although some success has been made in two dimensions for random surfaces. How- 
ever, when coupled to matter with central charge c > 1 these methods break down. 
Some observables like the intrinsic Hausdorff dimension are not accessible at all by 
analytic methods. In higher dimensions almost nothing is known analytically. The 
motivation of many of the conjectures and assumptions used for higher dimensional 
models is therefore derived from the experience obtained from the two-dimensional 
models. 

This work addresses these issues. We consider two models with different inter- 
nal local dimension, two and four. These models are namely hyper-cubic random 
surfaces, a discretization of bosonic strings, and four dimensional simplicial quan- 
tum gravity. To strengthen the common bases for the two models we discuss the 
subject in a somewhat broader context in the remaining part of this chapter. The 
motivation and definition of the models, considered in the numerical part of this 
work, is given in chapters |2| and |3|. 

1.1 Discrete random geometries 

Random geometries have attracted a lot of interest since it was realized that surfaces 
and higher dimensional manifolds are central to the understanding of the funda- 
mental interactions in nature. In this thesis we are interested mainly in quantum 
gravity and string theory, but they are also of interest in other fields like biology 

or condensed matter physics . All of these theories have in common that they 
describe statistical systems where we have to integrate over geometries. 

The thermodynamic properties of these models can be analyzed using standard 
methods of the theory of critical phenomena. It is generally believed that critical 
phenomena may play a similarly important role for random geometries as they do 
in statistical field theory. There it is known that close to a phase transition a 
certain degree of universality emerges in the sense that models with very different 
microscopic interaction exhibit the same critical behavior, described by a set of 
critical exponents. This is usually related to some common symmetry. Something 
similar happens for two dimensional quantum gravity coupled to matter with central 
charge c < 1. It was shown [Q] for this type of models that the entropy exponent 
7 depends only on the central charge but not on the type of the matter coupled to 
gravity^. In higher dimensions one numerically observes universal behavior in the 
sense that some phases of very different models have similar geometrical properties. 
For example both types of models considered in this thesis have a phase where the 
geometry is essentially described by a branched polymer. 

Discretization is a very useful tool in analyzing statistical field theories. It re- 
duces the number of degrees of freedom from non-countable to countable. A discrete 
lattice can serve as a regulator for UV-divergences and also opens the possibility for 
non perturbative approaches, such as numerical simulations. Typically, models of 
random geometries are invariant under coordinate transformations. Any regulation 
of the theory has to obey this invariance. Obviously the regulator introduced by 
discretizing on the geometric level is reparametrization invariant as this procedure 

^The definition of the exponent 7 is given in section [L3. 
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does not refer to any coordinate system. However, one has to ask the question what 
happens to difFeomorphism invariance, the basic symmetry of general relativity, in 
a discrete space. When taking the continuum limit of the discretized theory one 
has to check that this symmetry is restored, at least on large distances. 

Discretization can also help to define a measure. When performing the integra- 
tion over distinct physical states for the partition function the measure has to be 
chosen so as not to over-count surfaces by counting different parameterizations of 
the same geometry. For two dimensional quantum gravity it was shown that the 
critical behavior of the discrete dynamical triangulation model (a comprehensive 
introduction can be found in [||), where the integral over geometries involving the 
measure is replaced by a sum over triangulations, and the critical behavior of the 
continuum theory agree. 

The key ingredient that allows to go back from the discrete model to the con- 
tinuum is universality. If the discretized model is in the same universality class 
with the continuous version its critical properties describe continuum physics. The 
physical reason is that if the correlation length in units of lattice spacings diverges, 
for example at a second order phase transition, details of the discretization become 
unimportant and the system exhibits continuum behavior. 

1.2 Random walk 

A random walk describes the propagation of a point like particle in d - dimensional 
space. Despite its relative simplicity it contains many of the aspects common to 
other random geometries. As this model is well known also to readers from outside 
this field of research, it is perfect to demonstrate these aspects. For the details of 
the calculations we refer to 

In a quantum theory the amplitude G{x, x') for a particle to go from x to x' 
is given by the sum over all possible trajectories P^x' connecting these points. In 
Euclidean space we have 



where ^[Paja;'] is the classical action for the given path. As such we choose the 
simplest characterization of the path, its length: 




(1.1) 




(1.2) 



In the right-hand expression we introduced a parametrization 




(1.3) 



of the path. Expression ( |1.2| ) is invariant under reparametrizations given by 



(1.4) 
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with the function satisfying the conditions 

0(0) = 0, <A(1) = 1, ^>0. (1.5) 

We win show that with the action (|1.2| ) the quantum amphtude (^]^) describes the 
propagation of a free relativistic particle. Other terms in the action, for example 
higher derivative terms, can be relevant in the continuum as well. A non-trivial 
term is for example a coupling to the external curvature K. 

dl \K\ (1.6) 

For finite coupling e one finds the critical behavior of the ordinary random walk 
discussed below. In the limit e — > oo one finds a different class, smooth random 
walks [|,|. 

The expression ( |1.1| ) is rather formal. The measure T>P has to be defined in 
such a way that the integral counts only the number of distinct paths, i.e. paths 



which differ only by reparametrization (1^) have to be counted only once. In 



continuum formulations one formally defines this as the number of all paths modulo 
reparametrizations , 

VP 

A way to define the measure and at the same time regularize the model is to 
discretize it. This has to be done carefully, though: it is surely wrong to cut the 
interval < ^ < 1 into equally small pieces, as this would spoil reparametrization 
invariance. Instead of discretizing a certain parameterization one has to discretize 
the internal structure of the path. This amounts to considering only piecewise linear 
paths where each step on the path is of length a, i.e. we refrain from discussing 
structures smaller than a. This cut-off is invariant by definition. The possible length 
of the path is now I = na and the action is 5 = mo/. With this discretization the 
expression for the amplitude (^]^) is 

oo „ n 

Ga(x,x';mo) = ^exp(-moan) / de^J (^a ^ - (x - x')) , (1.8) 

n=l 4=0 

where the integral is over unit- vectors in W^. After Fourier transforming this and 
some algebra one finds 



with 



Ga{p;mo) = - (1.9) 

1 - exp(-moa)/(pa) 

f{pa) = j deexp(-iipe) ^ /(0)(1 - c^{paf + ■■■). (1.10) 

To recover the continuum from the discrete description the lattice spacing has 
to go to zero: a ^ 0. One gets the propagator for the free relativistic particle if at 
the same time the bare mass mo is renormalized such that 

/ \if \ 1 2 2 2- log(/(0)) 22 t^ ^^\ 

exp[—mQa)j[pa) — > 1 — c mpf^a , i.e. mo = he rUpf^a, (l-H) 
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where ruph is the physical mass. Putting the last three formulas together one gets 



Ga(p;mo) ^ ^ = ^r^Gcip;mph) (1.12) 

gZgz pz _|_ rn^f^ c^a^ 

when one collects the leading order a? terms. The discrete amplitude Ga in this 
limit is related to the continuum propagator Gc of the free relativistic particle by the 
divergent factor which can be interpreted as a wave function renormalization. 

It is interesting to note that these results are quite universal. They do not 
depend on the details of the discretization. For example one finds the same contin- 
uum limit for the random walk on a discrete lattice , where one has to sum over 
unit- vectors of in (LI) rather than integrating over the unit-sphere. One also 



finds universality concerning the action. Consider for example the covariant action 
of one dimensional gravity : 



S[x,g] = - I'd^yW) 
" Jo 



(1.13) 



where an independent internal metric "tensor" Qab^o, = 6 = 1 has been introduced. 
The partition function 

Z = l'^Vxe-''[x,g]. (1.14) 

contains in this case an integration over equivalence classes of metrics. The classical 
equations of motion obtained from ( 1.13| ) by variation of gab and x independently 



agree with those one gets from the action (|lj). But it is a non-trivial observation 
that they agree on the quantum level because in principle the metric field could 
induce different quantum anomalies. 



1.3 Random surface 

While point particles propagate along path a string sweeps out a surface when prop- 
agating. It is natural to expect that the generalization of the quantum amplitude 
( |1.1[) to surfaces describes the propagation of a free relativistic string if one chooses 
similar to (^]^) an action proportional to the area of the surface. Higher deriva- 
tive terms, like curvature-squared terms or extrinsic curvature terms may also be 
relevant in a possible continuum limit. It is, however, best to start with a simple 
action and check the relevance of other terms only afterwards. 
One such action is the Nambu-Goto action 



SNciSik)] 




(1.15) 

The metric of the surface is completely determined by the embedding of the surface 
x^ into a D-dimensional space. The 2-loop quantum amplitude is 

G{hM)= j VSeM-S[Si,,i,]), (1.16) 
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where the integration extends over ah surfaces Si^^^i^ which have li and I2 as a 
boundary. 

The non-polynomial action ( 1.15| ) makes the analytic treatment of the Nambu- 



Goto string very difficult. Therefore Polyakov suggested |^ a different string model. 
He introduced an internal metric Qab, a,b = 1,2 which is independent of the embed- 
ding. The action for this model is: 



S[g,x] =a d^C^/g 

JM{k) 



ah dx'^ dx^ 

9 + 



(1.17) 



The 2-loop Greens function now also integrates over the gab modulo reparametriza- 
tions : 

Giluh) = I ^ / VgXeM-S[g,x]). (1.18) 
The partition function for fixed area A is 

Z{A) = ['^1 V,xe-'^^^'^^5 (I d^Cv^ - A^ . (1.19) 

We define the entropy exponent 7 by the asymptotic behavior of the partition 
function for large A, 

Z{A) ^ e'^=^^^-^ (1.20) 

where fic is the critical coupling of the model. 

The classical equations of motion derived from the actions ( |1.15| ) and ( |1.17 ) 



agree because in the saddle point solution the metric gab takes the form induced by 
the embedding. But it is not at all obvious that the quantum theories are identical. 
It has been demonstrated [^] that the two quantizations are equivalent in the 
critical dimension D = 26. In lower dimensions it is not clear that the quantization 
'based' on the two actions are indeed the same. It is, for example, not excluded 
that the c = 1 barrier, which we discuss below, is absent for the model defined by 

In defining a string theory using the the Polyakov action one has to deal with 
the problem of avoiding over-counting of surfaces. The measure 

VfiO ^ ' 

is a functional integral over the space of metrics divided by the volume of the 
diffeomorphism group. In the continuum one can deal with this problem by gauge 
fixing. In the conformal gauge the metric gab can be written as 

gab = e''5ab (1.22) 

with only one parameter a. A naive argument which explains this result is that the 
four components of the metric tensor are not independent; a metric tensor is sym- 
metric, i.e. the off-diagonal elements are identical. Reparametrization invariance 
kills another two degrees of freedom, so one ends up with only one free parameter. 
This, however, is a bit too naive because the conformal gauge is accessible only if 
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the required coordinate transformation is non singular. It can be shown that on 
a sphere this is always the case. 

The action in ( |1.19| ) is invariant under the conformal deformation gabiO ~^ 
e'^^^^QabiC) of the metric. Therefore the classical action does not depend on a in the 
conformal gauge which reflects the triviality of the classical model. In the quantum 
theory, however, the measure T>gX transforms in a non-trivial way. This so-called 
conformal anomaly brings back the a dependence into the the effective action 

Si[(t)] = j d^z{2d,ad,a + ^?e''), (1.23) 



which is known as Liouville action. The gauge-fixing of the measure Vg introduces 

26 



a ghost 5; . The partition function for the Polyakov string in the conformal gauge 



IS 



Z = ypfjexpj-^^^ j d'^z{2d,ad,a + fi'^e'')Y (1.24) 



The in the Polyakov action ( 1.17 ) can be interpreted as D scalar fields cou- 
pled to the surface described by the metric gab- Note that with this interpretation 
the case D = 0, a surface embedded into a zero dimensional space, does make sense. 



In section |LJ we will show that this case describes pure gravity without matter. In 
it was shown that the entropy exponent 7 for two dimensional quantum gravity 
coupled to conformal matter with central charge c is given by 

^^,-l-V(l-,)(25_-c)^ (1.25) 

The central charge c in a sense counts the number of degrees of freedom of the 
conformal field theory. For a free bosonic strings it can be identified with the 
dimension of the embedding space. If one considers only unitary models with c < 1 
the exponent 7 takes only a discrete set of negative values 

1 

j = ,m = 2,3, .... (1.26) 

m 



Note that the formula ( 1.25| ) is not valid for c > 1 because it gives unphysical 



complex results; in the literature this is often referred to as the "c = 1 barrier". A 



more detailed discussion can be found in section 2.1 



1.4 Quantum gravity 

The classical description of gravity is given by Einsteins theory of general relativity 
(for an introduction see for example |l^, |l2[). The space-time metric describes the 
gravitational field. The field equations for the metric field are given by 

Rik — 7;9ikR = — j—Tik- (1-27) 
2 

where Gn = 6-6710^^^ m^kg^^s^^ is Newton's constant, gik the metric tensor and 
Tik the energy-momentum tensor of a matter field coupled to gravity. The Ricci 
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tensor Rif^ and the scalar curvature R = gikR^^ are functions of the metric tensor 
and contain first and second derivatives. The dynamical variable of the theory is 
the metric of space time. 

The classical equations of motion for this metric can be derived using the prin- 
ciple of least action, from the classical action. It is generally believed that it is 
possible to define a quantum theory of gravity using the path integral formalism. 
One assumes that not only the classical trajectory but all trajectories weighted by 
g-^la] contribute to quantum mechanical transition probabilities, giving 



S[9] 



(1.28) 



The integral is over all L'-dimensional Euclidean space time metrics which satisfy 
the boundary conditions g'^f\g^'^\ As usual in quantum field theory one has to 
switch to the Euclidean formulation. Unlike for conventional field theories it is 
still an open issue whether the Minkowski formulation can be recovered by a Wick 
rotation O]. In Z? < 4 dimensions the action in Euclidean space has the form 



S 



Gat 



\g\{2A-R)+S Matter, 



(1.29) 



where we have introduced the cosmological constant A. 

The quantum amplitude (1.28) describes a system of random geometries. We 
define the partition function of quantum gravity by 



Z 



Vf ' 



-S[9] 



(1.30) 



The case D = 1 has already been discussed in the context of random walks 
( |1.14| ).The distance I = J y^j^fdC is an invariant internal characteristic of the met- 
ric, but we cannot define an internal curvature. Therefore the scalar curvature R 
term in ( 1.28] ) drops out. Comparing with the action (1.13) we see that the two 
expressions differ by the term 



S[x,g] 



1 



a 



dev5(0 



(1.31) 



which describes the embedding of the random walk in a A-dimensional space. For 
gravity is does not make sense to interpret the x'^, /x = 1, . . . , A as an embedding, 
but one can interpret this term as the contribution of A scalar fields. The case 
A = 0, i.e. the random walk embedded in a zero dimensional space, describes one 
dimensional quantum gravity without matter fields. 

For D = 2 quantum gravity is described by the Polyakov string ( 1.1^ ) embedded 
in zero dimensions. Note that for surfaces the curvature i? is a local invariant, so 
we have a new degree of freedom compared to the one dimensional model. On the 
other hand it is known from the Gauss-Bonnet theorem. 



d^^ 



\g\R = Attxe, 



(1.32) 
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that the integrated curvature is a topological invariant with xe = 2 — 2h, where 
h is the number of handles in the surface. The Euler character is xe = 2 for 
the sphere, xe = for a torus and so on. For fixed topology the curvature term 
m the action ( pTloD gives therefore no dynamical contribution, which makes the 
classical model trivial. One might suspect that this is also the case for the quantum 
model. In fact the two dimensional model does not contain gravitons, however the 
conformal quantum-anomaly generates the non-trivial behavior discussed in the 
previous section. 

Quantum gravity in physical four dimensions is one of the models considered 
in the numerical part of this work. Therefore we discuss this case in detail in 
chapter ^. Here we only want to mention that for the time being no continuum 
solution for the model is known. The discrete approach, which is very successful 
in two dimensions can be generalized to four dimensions. With this method at 
hand numerical simulations allow us to survey properties of quantum gravity in 
this interesting case. 



1.5 Matrix models 



Figure 1.1: A part of a M^-Graph. The dotted line depicts the triangulation dual 
to the graph. 



In dynamical triangulations [14, 15, [T^, the path integral ( 1.3C| ) over two dr 



mensional equivalence classes of metrics is replaced by a sum over triangulations 



/ 



The symmetry factor cr(T) is a discrete remnant of the reparametrization group. 
It is mainly important for small triangulations because typically larger triangula- 
tions are not symmetric. A triangulation is a surface obtained by gluing together 
equilateral triangles along their edges so that exactly two triangles meet on each 
edge with possible exceptions on the boundary. As in the continuum we choose the 
action to be proportional to the area. In the discrete model the area is given by 
the number n of triangles in the triangulation. 
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In the partition function 



Z{^l,G) = i^exp(-/xn) = ^exp | J]exp(-/xn)A4^)(n). (1.34) 

T X " 

the topology is not a priori fixed. However the number of graphs with free topology 
grows factorially with n and therefore the naive topological expansion in the number 
h of holes in the surface is as stands ill defined. Therefore calculations and numerical 
simulations are in most cases done for ensembles with fixed topology. At the end 
of this section we discuss the double scaling limit, which can be used to define the 
model with coherent contributions from all topologies. 

The dynamical triangulations model can be formulated as a matrix model, de- 
fined by the generating function 

Z = log j dM exp{-TrM'^ + -j=TrM^), (1.35) 

where M is a x hermitian matrix and dM the fiat measure 

N 

dM = J\dMii d{ReMij)d{ImMij). (1.36) 

i=l l<i<j<N 

Doing perturbative expansions in g one gets Feynman diagrams similar to those 
depicted in figure |1.1| . Note that the dual graph, drawn in dotted lines, is a tri- 



angulation. Because ( 1.35| ) generates all possible graphs it sums over all triangu- 



lations. It should be emphasized that the results presented below do not rely on 



using triangulations. In |17] it was shown that the continuum limit for the model 
with a lattice built from squares, pentagons, etc. agrees with the result obtained 
for triangulated surfaces. 

Differently from a standard Feynman graph, the propagator 

i ^ 1 

< MijM^i >= 6u6jk = . 

J ^ k 

carries two indices. This is necessary to interpret the graph as a two dimensional 
surface; a priory Feynman graphs do not have any "dimension". Moreover, ( |1.35D 
generates orientable surfaces because M is hermitian. Due to the Wick theorem 
one gets the nth term of the perturbative expansion by summing over all diagrams 
with n vertices, where each vertex has a weight 




(1.37) 



A closed index loop contributes a factor A^ because the trace sums over the indices 
of the matrix. A diagram D with V = n vertices, L closed index-loops and P 
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propagators therefore contributes 



^toiV-V^,-'iV-=,-'_A.«. (1.38) 

where a is the symmetry factor of the diagram. The expression on the right hand 
side is obtained by noting 'iV = 2P and using the Euler theorem V + P — L = 
2 — 2h = Xi where x is the Euler character of the surface, which uniquely identifies 
the topology of the graph. 

In the diagrammatic expansion the partition function of the one matrix model 
reads 

^{g,N) = Y,mY.9'' E ^ = E^"E5"4^^N- (1-39) 

h n D&Vh^„ ^ ' h n 

The number J\fx\n) of graphs can be calculated with the method of orthonormal 
polynomials one finds 

A4'^(n) ocn^^-^^-". (1.40) 

The entropy exponent 

7x-2 = cx, c = -^ (1.41) 

is linear in x- The constant Qc-, the critical coupling, does not depend on the 
topology. 

In the limit N ^ oo, the so-called planar limit, only diagrams with spherical 
topology X = 2 contribute to the sum over topologies in the partition function 
( 1.39| ). Contributions from higher genus diagrams are suppressed by N~^^. On the 
other hand these contributions are enhanced as g ^ Qc, which can be seen from the 
asymptotic behavior of the partition function with fixed topology : 

Zxia) = Y,g^M';^\n) ^{g- g.f-'-^ (1.42) 

n 

This suggests that if the limits N ^ oo and g ^ gc are taken together in a correlated 
way the large genus suppression can be compensated by the g ^ gc enhancement. 
This would result in a coherent contribution from all genus surfaces p|, 
This so-called double scaling limit is defined by taking the limits N ^ oo, g ^ gc 
while holding fixed the the 'renormalized' string coupling 

= N{g - g,)-~^. (1.43) 



By comparison of (1.39) with the partition function of dynamical triangulations 
(|1.34|) we find the relations 

iV = exp ^, g = exp{-n) (1.44) 

for the coupling constants. For spherical topologies we can compare the exponent 
72 with the the continuum results (1.25). The value 72 = — ^ obtained from (1.41) 
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Figure 1.2: A part of a hyper-cubic surface and the representation of the surface 
as a branched polymer 



agrees with the value 7 = — ^ obtained in the continuum for pure gravity, i.e. c = 0. 
Some variants of the matrix model discussed here can be used to describe discrete 
gravity coupled to matter. For example a particular two matrix model describes 
the Ising model (c = ^) coupled to gravity |22]. One finds 7 = — | in this case 
which is in agreement with the results obtained in the continuum. 



1.6 Branched Polymers 

Many models of random geometries, for example the two types of models considered 
in this thesis, have a phase which can essentially be described by branched polymers. 
Following Q we want to derive the result 7 = 1/2 for the entropy exponent because 
we frequently refer to this result. As an example for a branched polymer structure 



we show in figure 1.2 a small hyper-cubic surface, which can be decomposed into 
small components by cutting open loops of length two. The individual components 
are not critical, their size is of the order of the lattice spacing. The only way 
the number of plaquettes can grow to infinity is by successive gluing of branches, 
all the dynamics lies in this gluing. This dynamics is described by the branched 
polymer model : in the right hand side of figure |1.2| the individual components are 
represented as links li with an associated chemical potential fi. We also introduce 
a weight fn for the joining of n links li at a vertex Vn- The partition function for 
the branched polymer model is 

BP 

where the sum runs over all branched polymers, i.e. tree graphs if we restrict 
ourselves to spherical topology. To extract 7 we consider the one-point function 

G{fi)^c{fi-ti,)'-\ (1.46) 

the amplitude of rooted trees with one marked vertex. Such a tree can be built 



recursively as seen in figure 1.3, where the blobs represent a rooted tree. This can 



be used to write down the self-consistent equation [23| 



= e-f^ (1 + f2G^ + hGl + ...). (1.47) 

Solving for gives : 
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Figure 1.3: The self-consistent equation for rooted trees. 



J2 ^, ji ^_ = F{G^) (1.48) 



G^ 



We conclude that the lowest value /ic of fi for which eq. ( 1.48|) has a solution is the 
minimum of F[G^) if the weights /„ are positive : 

F{G,)\,=,^=Q (1.49) 



dG^ 



Therefore the linear term drops out if one expands (/i — ^c) in G^ : 

fjL- Hc^ [Gf,^ - G^f i.e. G^ w G^^ - c^/ ii - fj-c- (1.50) 
Comparing with eq. ( |1.46D one finds 7 = ^• 



1.7 Monte Carlo methods 

The algorithms used to simulate the two different models, hyper-cubic random sur- 
faces and four dimensional simplicial quantum gravity, are both based on the same 
numerical technique, namely Monte Carlo methods. We give a detailed description 



of the algorithms separately for each model in section |2.5| and |3.5| . Here we want 
to summarize briefly some well-known facts about Monte Carlo algorithms which 
we use as a foundation for the more detailed discussion. 

Despite the apparent simplicity of a discrete partition function in the form 

Z = ^exp(-5H) (1.51) 
r 

it is only for a few simple systems possible to actually perform the summation. 
Even with a computer it is only for very small systems possible to sum over all 
states a G T because the number of states MaiV) in general grows exponentially 
with the volume. However, only comparatively few states give a significant contri- 
bution to the partition function. The overwhelming majority of the a are either 
entropically suppressed or damped by the Boltzmann factor exp{—S[a]). Monte 
Carlo algorithms are a method to sample the important contributions by choosing 
at random n samples cTj G T with the probability distribution 

7r((Tj) oc exp — S'ftTj]. (1.52) 
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This is the thermodynamical equilibrium distribution. One can therefore approxi- 
mate observables by their expectation value 



< O >= (1.53) 

n 

It follows from the central limit theorem that the error depends on the number n 
of measurements as cr{0) oc l/\/n if the cjj are independent. 

A problem is of course how to choose the cTj at random from the state space 
r. Dynamical Monte Carlo methods simulate a Markov process which generates a 
sequence of ai. The main ingredient of the Markov chain is a set of transformation 

C7i (Tj+l (1.54) 

and a probability matrix iJ^{(Ta-,(^b))ab- One can show that the n-step transition 
probability p(")(cr ajO'ii) converges to a equilibrium distribution 

lim ^ 7r(crb) (1.55) 

t—^00 

if the following conditions are satisfied: 

• Ergodicity For each pair <ti,(T2 G T exists a sequence of transformations 
which as a whole transform ai into (72. 

• Stationarity The following equation holds for the components of the tran- 
sition matrix and for all ah GT: 

Vr((Ta)7^(<Ta, CTfe) = 7r((Tfc) (1.56) 



The solutions of the detailed balance equation 

Tr{aa)V{cra,(Jb) = iT{ab)V{crb,(Ja) (1.57) 

satisfy equation ( 1.56| ). This condition is stronger than ( 1.56 ) but practically more 
convenient. In praxis the following solutions of this equation are important. With 
the heat bath weight 

V{aa,c7b) = 7r{ab) (1.58) 

the transition probability to the new state ah is independent from the old state 
The Metropolis weight accepts all transformation for which 7r(6) > 7r(a). The 
weight for the inverse transformation can be deduced from ( |1.57D . The shorthand 
notation for both directions 

P(f7„afe) = mini 1,4^1 • (1-59) 
is especially useful for a computer implementation. 
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Chapter 2 

Hyper-cubic random surfaces 



2.1 Introduction 



The Nambu-Goto action ( [L.15 ) describes a bosonic string embedded in a di- 



mensional space TZ. The world-sheet metric is, different from a theory of quantum 
gravity using the Polyakov string, completely determined by the embedding of the 
string in TZ. One can therefore attempt to discretize the model by discretizing the 



embedding space. This idea was used long ago in \24\ to define a model of hyper- 
cubic random surfaces. However, in [^] it was shown that this discretized model is 
trivial because the surfaces which dominate the partition function have the struc- 
ture of branched polymers, basically one dimensional highly branched thin tubes. 
It was generally believed on grounds of renormalization group arguments that one 
cannot cure this behavior by adding local terms to the action. It therefore came 
as a surprise when some years later numerical results for a slightly modified model 
with a local constraint indicated non-trivial behavior |2^. In D = 4 the anomalous 
scaling dimension « 1 was measured to be different from zero and for the entropy 
exponent 7 ~ 0.25 was obtained. 

In today's language of random surface theory one can interpret this model as 
a random surface model coupled to matter with central charge c = D > 1. The 



KPZ-formula ( |1.25 ) classifies 2(i-gravity coupled to matter with c < 1. The c > 1 



regime however, where ( 1.25| ) gives an unphysical imaginary exponent 7, is still not 



completely understood. Discretized models of 2d-gravity are, on the other hand, 
well defined for c > 1. Numerical simulations clearly indicate branched polymer 
behavior for c ^ 5, although in the regime 1 < c ^ 5 the numerical evidence is not 
so obvious. They seem to indicate a smooth cross-over to the branched polymer 
phase as c increases. Based on the numerical evidence one might conclude that 



there is a different critical behavior for 1 < c ^ 5, although it has been argued |2J], 
based on a renormalization group analysis, that the universal behavior for c > 1 
models is always a branched polymer structure. This behavior, however, is masked 
by the presence of the c = 1 fixed point, which becomes complex for c > 1. This 
would lead to exponentially enhanced finite size effects. A recent numerical work 
[p7| supports this conjecture. 

In view of this result 7 ~ 0.25 for the modified hyper-cubic random surface 
model, considered in [p^] is, if true, very interesting ||2^. It would provide a very 
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simple model with a non-trivial value, besides the mean field result 7 = |, of the 
series 

7 = -,n = 2,--- (2.1) 

n 

of positive 7 values discussed by Durhuus However one may ask how a local 

constraint, which disappears after one coarse graining step in a renormalization 
group analysis, could influence the large scale properties of the model. On the 
other hand, it is just the fact that the coarse grained surfaces in general does not 
belong to the constrained class of surfaces which invalidates the triviality proof 
[ p5[ for the unrestricted model. In addition a different universal behavior would 
indicate the breakdown of universality because in a sense the constraint is nothing 
but a detail of the discretization. 

A mechanism which might induce non-trivial behavior into the model is sug- 
gested by the renormalization argument used to derive the series ( |2.lD . In this 
argument the surface is decomposed into shortest loop irreducible components, i.e. 
components which cannot be cut into two parts along a (model dependent) shortest 
possible loop. For systems with n > 2 the entropy of the surface is still dominated 
by branching. But differently from branched polymers the components of the sur- 
face which form the branches are themselves critical with 7 < 0. Such a behavior 
is induced by matter field(s) with the couplings tuned to their critical values. The 
local constraint introduced into the model (2.5) can be written as an additional 
term in the action (|2.6| ) as a coupling to a external curvature. For the (actually 
infinite) coupling, introduced by the constraint, the model could be critical and 
above scenario could apply. 

We will address the question if a local constraint induces really a different critical 
behavior. We furthermore generalize the unrestricted hyper-cubic random surface 
model by adding to the action a term coupled to the external curvature of the 
surface. This model allows an interpolation between the |25| unconstrained and 
the potentially non-trivial [26| constrained version of the model. In a optimistic 
scenario one may hope that by tuning the couplings one may induce a non-trivial 
phase structure with a weak coupling branched polymer phase and a strong coupling 
fiat phase. At the phase transition the model might have a another value of 7, just 
as in the scenario described by Durhuus. 

^see also section ^.4.] 
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2.2 The model 



A hyper-cubic surface is an orientable surface embedded in Z^, obtained by 
gluing pairwise together plaquettes Px^y along their hnks until no free link is left. A 
plaquette Px^y is a unit-square which occupies 
one of the unit-squares , /x = 1 • • • in the 
embedding lattice. In the following we use the 
word square to refer to the embedding lattice 
and the words link or plaquette to describe the 
internal connectivity of the surface. A link is 
shared by exactly two plaquettes, which are 

glued together along each one of their edges. Figure 2.1:' A part of a hyper-cubic sur- 
Note that a link or a square in the embedding face in D=3 dimensions 
lattice can be occupied more than once. In other words, the surface is not self 
avoiding; touching and self intersection are allowed. 

The hyper-cubic random surface model is defined by the partition function: 




2(/3) 



E&S 



(2.2) 



where the sum runs over an ensemble S of hyper-cubic surfaces with spherical 
topology. In the large volume limit, the number N{A) of surfaces with a given area 
A is expected to grow like 

M[A) ^ e'^=^A^-^ (2.3) 

where 7 is the entropy exponent. Let us for the moment consider the action with 
only one term: 

Sif,) = M. (2.4) 



The partition function ( |2.2D is well defined only if the coupling /i is larger than the 
critical value /ic- The continuum limit is obtained by taking the lattice spacing a 
to zero and at the same time fi ^ fi'^. 

There are three possible ways to embed two neighbor plaquettes pi,P2 in a 
hyper-cubic lattice Z^. The external angle @i, which we assign to link /, can 
therefore take three possible values: 



0: 



TT . 

2 ■ 



tt: 



pi and p2 occupy two neighboring squares in the same coordinate plane, 

pi and p2 occupy two neighbor squares in different coordinate planes, 

pi and p2 occupy the same square. We call such a configuration self-bending. 



The model defined so far describes what we called the unconstrained model. 
The constrained model considers the ensemble of hyper-cubic surfaces 



5' = {5 e ^l/S does not contain self bendings}. 
This restriction can be formulated by introducing a potential 



(2.5) 



(2.6) 
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to the action with an infinite coupling e. The angle @i is the angle of two plaquettes, 
which have the edge I in common, measured in the embedding space. Therefore 
this term couples to the external curvature of the surface. 

It seems a bit artificial that this potential only contains contributions from links 
I with @i = TT. We therefore also use the potential 

si^tiEM ='-^il-cos&i) = eY^ (s^,e, + %2,e) • (2.7) 

In summary: The action of our model is 

5[E](/x,e)=M + 5«[^], z = 0...2 (2.8) 

where 

^S(e) = (2.9) 

formally denotes the action without a curvature term. The canonical partition 
function is 

Z(A6)= E e^^I^K^'^), (2.10) 

where the action S contains one of the external curvature terms discussed above. 
The sum runs over hyper-cubic- surfaces with fixed area A. 
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2.3 Universality 



In this section answer the question of whether the constrained hyper-cubic random 
surface model has a critical behavior different from the unrestricted model as in- 
dicated by earlier numerical results [^6| in four embedding dimensions. With the 
algorithms and the improved computer hardware available today we can simulate 
systems up to two orders of magnitude larger then those used in earlier works. 

The model is defined by the partition function (|2.2| ) where the summation runs 
over the constrained ensemble (^]^). The action ( |2.8| , i = 0) contains no external 
curvature term. We measure the universal entropy exponent 7 defined in (|2.3| ), 
which is a good indicator for different phases of the model. For branched polymers 
it takes the value 7 = ^- A different value indicates a different geometry. This 
will become clearer in section p. 3. 2 where we use a geometrical decomposition to 
measure the exponent 7. 

2.3.1 Old method 

As the first step we want to make sure that our program, which uses an algorithm 
different from the one used in the earlier works, gives the same answers if one uses 
a similar numerical setup. We simulate the grand canonical ensemble to extract the 
distribution Nf^{A) of surface areas measured with coupling fi. Taking into account 
the Boltzmann factor one can fit N^{A) to the theoretical distribution M{A) eq. 
( [2.31 ). Besides an uninteresting normalization constant one gets an estimate for 7 
and t = fic — fJ-- However, one has to fine-tune fi very precisely so that t = 0, because 
for large {tA) the distribution is dominated by the leading exponential factor and 
therefore the sub-leading behavior is difficult to observe in this case. We have used 
this method anyhow for small surfaces and got a value for 7 which is compatible 
with the results presented below. 

In |2^] a clever way to avoid this problem is used. The unknown parameter fj,c 
can be absorbed into a constant if one divides two histograms N{Ai),N{A2) with 
6 = Ai — A2 kept constant : 

^1 _^/.2g(M2-m)Ai f £1 \ ^ (2.11) 



^2 

In this case the factor 



i2 



gM2 ^ g(Mc-M2)<5 (2.12) 

is just a constant. 

In the numerical simulations we use our Monte Carlo program to generate the 
distributions N^(A) for fi < fic- We prevent the surface from growing to infinity 
by introducing an upper bound A^ax for the area. This does not spoil the detailed 
balance condition but ergodicity may in principle be affected for surfaces not much 
smaller then Amax- In the verification of the algorithm (section p. 5.1 ) we, however, 



did not observe any effect on the distribution of surfaces induced by the upper 



bound. Following [26| we have anyhow discarded the data for surfaces of size 



[-Amax — 20, Amax] when analyzing the data. We have measured Nfj_{A) with two 



different Amax and for different couplings. In table 2.1 we summarize the number 
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■^max — 300 


Amax = 5000 


fj- = 


1.60 




3 X 10"^ 


^J' = 


1.65 


1 X 10^° 


3 X 10^° 


M = 


1.70 


2 X 10^° 


3 X 10^° 


^ = 


1.75 


1 X IQio 





Table 2.1: Monte Carlo statistics in iterations collected for parameters Amax and 



P(Y) 

0.20 - 



0.10 - 



0.10 0.20 0.30 0.40 

Y 

Figure 2.2: The distribution of 7 values 



of iterations collected for and /U. After dividing each 

dataset into three parts there are 9*9 = 81 ways to combine data; for -A^T^ax — 5000 

We have extracted 



in formula (2.11) with data for A. 
combinations. 



300. 



7 



for the 81 



In figure 2.2 we show the distribution of 7 values obtained with 
this procedure. Taking the average we get 7 = 0.26(2). This is in agreement with 
the result 7 = 0.24(3) in |2£]. We have repeated the procedure for the model 
which allows for self-bending surfaces and found 7 = 0.44(3) which is in reasonable 



agreement with the expected result 7 



2.3.2 The Baby universe method 

The baby universe method [^] has become a standard tool for measuring the 
entropy exponent 7 in numerical simulations of random surfaces. It is very conve- 
nient because it does not require fine tuning of parameters and allows the use of 
large systems in the simulation. In this section we present the numerical results 
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A-B 



B 



Neck 



obtained with this method for hyper-cubic random sur- 
faces with up to 80000 plaquettes. 
A baby-universe is a part of a surface S with the topol- 
ogy of a disc. It is connected with the rest of S in such 
a way that the disc-perimeter / is much smaher than 
the disc area B: f < B. In figure BXl this is depicted 



Figure 2.3: A baby surface 



graphicahy: a baby universe is a part of S connected to 
the rest of 5 by a bottle-neck. If one looks at this neck 
from far enough away one sees it as a point at which two surfaces are joined, one 
with area B, the other with area A — B, where A is the area of the whole surface 
5. 

One can approximate the baby universe distribution |33] by : 



NAiB)^Mi{B)MiiA-B) 



(2.13) 



where Mi{A) is the number of discs with loop / and area A. This number can 
be expressed by the number of spherical surfaces M{A) eq. ( |2.3| ). Consider for 
example the disc with area A and a unit-square as the boundary. Such a disc can 
be created by removing one plaquette from a spherical surface with area A + 1. 
Taking into account the factor (^ -|- 1) which counts the number of plaquettes that 
can be removed from the sphere one gets 



MiiA) = (A + 1)AA(A + 1). 



(2.14) 



This formula is approximately true also for different boundary shapes. Substituting 
this expression and ( p. 3D in ( |2.13| ) gives 

Na{B) ~ c{A){B{A - B)y-\l + corr.) (2.15) 

for the distribution of baby surfaces with area B on a surface with total area A. The 
leading exponential contribution c{A) = cexp {fJ-c{B + {A — B})) depends only on 
the area A of the whole surface and not on B. 

This is very convenient in numerical simulations where one can use ( p. 15 ) to 
measure the exponent 7, which is part of a volume correction term in the canonical 
(fixed area) ensemble. This means that one does not need to fine-tune fi to /ic- 
However, for technical reasons, we still have to tune with comparatively low pre- 
cision the coupling fi to its pseudo-critical value fipc- Our algorithm is inherently 
grand canonical since ergodicity requires volume fluctuations. We extract the infor- 
mation for the canonical ensemble with A = Aq from the simulated grand canonical 
ensemble by measuring only if ^ = ^q- Inefficient large volume fluctuations about 
Aq are suppressed by a potential 

U{A) = ^-{A-Aof (2.16) 

which we add to the action (|2.8| ) with an appropriate choice of 6. We used 6 = 0.002 
in the simulations. For the simulation the coupling n has to be tuned so that 
< A >~ Aq. The Gaussian form of the grand-canonical ensemble allows us to 
estimate /i — fipc by measuring < A >, i.e. 



H - UpMo) « 5 (< ^ > -Ao) 



(2.17) 
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Figure 2.4: Pseudo critical coupling Hpc{A) 



We use this formula to improve, recursively, the estimate for /i in the simulations. 



In figure 2.4 we show the pseudo-critical coupling obtained with this procedure. 



We choose the coordinates 

y = ln[iBiA-B)fNAiB)] 

X = ln[B{A-B)] ^ ' 

to analyze the measured baby universe distribution Na{B). In these coordinates 
the exponent 7 can be extracted from the numerical data with a linear fit to the 



equation y = + b. In figure 2^ we show a typical baby universe distribution 
using the coordinates (2.18). The total surface area is ^ = 20000. For small values 
of B we see clear evidence of finite size effects; the data-points do not lie on a 
straight line. On the other hand for large values of B one can see strong statistical 
noise which comes from the low frequency of large baby universes. For example, for 
A = 10000 the probability of finding a baby universes of area B = 4999 is 7 x 10~^. 
When estimating the value of 7 one has to deal with the finite size effects. In our 
analysis we discard data for small baby universes. To decide, which part should be 
kept, we look at the behavior of jeffiB) obtained from the fit to the data for baby 
surfaces larger than i? as a function of B. As one can see in figure (p.6|), where 
we show as an example 7e//(-B) measured for A = 20000, that jeff{B) grows to a 
plateau value. We use this value to estimate 7e//- 

We analyzed the baby universe distribution for different canonical volumes. In 
figure (2^) we show 7e// as a function of the canonical volume A. For comparison 



we also performed simulations for the unrestricted ensemble. As discussed above, 
for this case one can show analytically that 7 = ^- For the constrained model 
without self-bendings (lower curve) one observes stronger finite size effects. The 
asymptotic value of 7e//, however, is the same. We conclude that the constrained 
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Figure 2.7: The effective entropy exponent "jeff as a function of the canonical 
volume A for ensembles including (upper curve) or excluding (lower curve) self- 
bending, respectively. 



hyper-cubic random surface model possesses the same critical behavior as the un- 
constrained model. This means that the model is universal in the sense that the 
details of the discretization are unimportant. However the result 7 = ^ contradicts 
the results obtained with the method used in [26|, presented in the previous section. 
The reason for this contradiction is easy to understand if one looks at figure 2.7 . 
In the last section surfaces with an area 200 < A < 280 were used to extract an 
estimate for 7. From figure 2.7 one sees that for surfaces of this size the effective 
7gjj PS 0.25; on such small surfaces the finite size effects are still dominating. 



2.3.3 Discussion 

We have measured the value of the string susceptibility exponent 7 for hyper-cubic 
surfaces embedded in a 4 dimensional target space for an ensemble of spherical 
surfaces with and without self-b endings. In both cases we find numerically, for 
large enough surfaces, 7 ~ 1/2, which is compatible with a branched polymer 
behavior. The local constraint first discussed in |^6| is not sufficient to change the 
critical behavior of the model. However, in the next section we will show that a 
generalized external-curvature term can change this behavior. 

It is intriguing to note that the discretized model of Nambu-Goto strings suffers 
from the same pathological behavior as that of Polyakov string theory. It collapses 
down to a branched polymer, i.e. essentially a one-dimensional object. This result 
suggests that both models have the same universal behavior. Unfortunately the 
Nambu-Goto string is not defined in D < 2, so it is not possible to compare the 
discretized string models in D < 1 where one finds non-trivial behavior for the 
Polyakov string. 
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2.4 Hyper-cubic random surfaces with extrinsic curva- 
ture 



2.4.1 Introduction 

The branched polymer behavior of the constrained hyper-cubic random surface does 
not rule out the possibility that one finds non-trivial behavior for the model with 



the generalized extrinsic curvature action (2/7). With infinite coupling e the generic 



surface is, up to finite size effects, flat. For finite couplings one may hope to find a 
transition from the branched polymer to a flat phase. 

However it was shown in [^3| for the model with extrinsic curvature that "the 
critical exponents take their mean field value if the susceptibility of the model and 
a coarse grained version of it both diverge" . In other words under the assumptions 

7 > and 7 > 0. (2.19) 

the model is always in the branched polymer phase. This result was obtained with 
a renormalization group argument where the branched polymer was decomposed 
into "blobs" , components which can not be cut into two parts along a loop of length 
two. The exponent 7 is the entropy exponent for these blobs. The equivalence of 
( |2.19| ) with the statement given above can be seen from the expected behavior of 
the susceptibility 

x(m) ^{^^- f^c)-\ (2.20) 

which diverges for fj, ^ Hc only if 7 > 0. 

A key step in the derivation of this result in |p3i is the formula 



x(/i) = (2.21) 

which relates the susceptibility x(/^) of the original model to the susceptibility xifi) 
of the decomposed blobs with renormalized coupling fi. Obviously x = 1 if the 
susceptibility x diverges. Under the assumptions stated above, i.e. 7 > for the 
blobs one concludes that they are not critical p, > fic and x is analytic at this point. 
With a Taylor expansion one gets the self consistency relation 7 = 1 — 7 which is 
solved by 7 = ^, the generic value for branched polymers. 



However, Durhuus emphasized in |30] that the condition x = 1 alone does not 



imply that the blobs are non-critical. If 7 < the susceptibility x does not diverge 
and X = 1 can be satisfied at the critical point fl = flc of the coarse grained system. 
At such a point the entropy is still dominated by the branching of the surface but 
the branches themselves are critical, which in effect changes the exponent 7 for the 
whole system. If one assumes that 7 takes the KPZ-values ( |1.25| ) one can derive 
from (|]2l|) the series yU.u). 

The derivation of that formula is rather formal. It does not provide a description 
how to construct a system which has this property. We want to check numerically if 
the hyper-cubic random surface model exhibits non trivial behavior at an eventual 
crumpling transition. In other words we check if the assumptions ( |2.19| ) are fulfilled 
in the whole range of couplings e. Non-trivial behavior can be expected only if these 
assumptions are not satisfied, which we expect if e = ec is at a critical point. 
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Figure 2.8: The susceptibility Xext in D = 3 dimensions 
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Figure 2.9: The susceptibility Xext in D = 4 dimensions 
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Figure 2.10: The distribution of the average external curvature in D = 3 dimen- 
sions obtained with re-weighting. The non zero probabihty for unphysical negative 
values < is an artifact of the extrapolation to this region in the re-weighting 
process. 



2.4.2 Numerical results 

The simulations are done in the same quasi-canonical way already discussed for 
the constrained model eq. ( 2.16| ). To localize possible transitions we measure the 



energy fluctuation of the external curvature field 

Xe.t = log Z{A, e) = \ (2-22) 

and search for a peak in this observable. The average external curvature is propor- 
tional to the first derivative of the free energy : 

< ® >= ^1 ^) = 2^ 5: E 0^^-''"^- (2-23) 

s I 

We simulated in both three and four embedding dimensions and scanned the 
coupling range < e < 3 for peaks of Xext- In both cases we found a single peak. 



We used re- weighting methods |34] to extract the shape of the peaks shown in 
figures |2^ and |2]^ from four independent measurements at e epc per volume. 

The peaks grow quickly, in fact a bit faster than linearly, with the volume. 
This indicates a first order phase transition, which is confirmed by a look at the 
distribution of external curvature 0, where the bar indicates the average taken 



over a given lattice. In figure 2.10 we show this distribution for D = 3 and the 
coupling e epc close to the pseudo-critical coupling. One observes a clear signal of 
a first order phase transition, namely two maxima separated by a minimum which 
becomes deeper as the size of the surfaces is increased. This clearly indicates two 
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Figure 2.11: The pseudo-critical coupling ecrit{'^/A). 



separate phases. In one phase the average external curvature is close to zero, the 
typical surface is flat. 

To estimate the critical coupling we use that for a first order phase transition one 
can expect 

(2.24) 



(^) = + ^ + o( j,; 



for the scaling of the pseudo-critical couplings €pc- In figure |2.1l| we show the 



numerical estimates for e^c- Note that the pseudo critical coupling decreases with 
the volume which indicates that ec is finite in the thermodynamical limit. With a 



linear fit to equation (2.24) we find for the critical coupling 



D 
D 



3 tc 

4 Er. 



1.530(3) 
1.643(4) 



To demonstrate the change in the geometrical behavior we consider the radius 
of gyration 

G\A, e) =< i 5](5^(p,,,) - S,{po)f >, (2.25) 



A 

Px.y 

which is the average squared distance of the plaquettes to some reference plaquette 
Pq. The are the coordinates of the plquette in the embedding space. In the 
numerical simulations we do a stochastic average over the pq by choosing about 10 
% of the plaquettes at random as po per measurement. One can use the radius of 
gyration to define the external Hausdorff dimension dh 

G^{A,€) ^ A^/'^'^. (2.26) 

In a sense the Hausdorff dimension is the largest possible dimension for the embed- 
ding space which can still be completely filled by the surface. One has for example 
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Figure 2.12: The scaling of the radius of gyration. 
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Figure 2.13: The phase diagram for the Hyper-cubic random surface model with 
extrinsic curvature. One finds qualitatively the same behavior in D = 3 and D = 4 
embedding dimensions. 
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G oc A for a flat surface and therefore dh = 2. In figure 2.12 we show G{A,e) 
measured in the branched polymer phase, at e = 0, and in the strong coupling 
phase, at e = 2.0. With a fit to (|]2|) we extract the Hausdorff dimension in the 



two phases: 

e = 0. dh = 4.07(3) 
e = 2. dh = 2.00(3). 

The value found for e = is two standard deviations away from the expected result 
76| but this should presumably be attributed to finite size effects. 



The results are summarized in the phase diagram figure |2.13. The model is 



defined only for fi > /Xc(e) which defines the critical line depicted by the solid line 
in the phase diagram. Above this line we find two phases. For e < Ccrit the system 
is in the branched polymer phase with d^ ~ 4. If we assume hyper scaling 

u = d],\ 7 = z.(2-7?) (2.27) 

the anomalous scaling dimension is = as expected for a branched polymer. The 
strong coupling phase is flat with dh = 2. In this phase basically no baby universes 
are present. For example, on a surface of size A = 500 we found at e = 2 on average 
only 0.06 baby universes of size B = 5. The largest one found in this phase is of 
size 33 and appeared only once in 200000 measurements. For comparison in the 
branched polymer phase at e = we found in average 10.6 baby universes of size 5 
on the surface, the probability for the largest possible one with B = 249 is about 
0.02. Therefore one can not use the baby universe distribution to measure 7 in the 
flat phase. However the fact that there effectively are no baby surfaces means that 
7 is large negative. Therefore the existence of this phase is not in contradiction 
with [^] because the assumptions ( p. 19 ) are not met. 



We hoped to find 7 < for the branches at the critical point. Instead we found 7 
for the whole system large negative in the strong coupling phase. At the transition 
point the exponent 7 is not well defined because the transition is first order with 
different values for 7 in both phases. The close to horizontal critical line in the 



phase diagram fig. |2.13| for the flat phase reflects the fact that the curvature field 
is almost frozen at < >= 0. A increase of the coupling e has no further effect on 
the geometry and therefore also no effect on the entropy log(/ic) of the surface. 

2.4.3 Discussion 

We have investigated the phase diagram of hyper-cubic random surfaces with an 
extrinsic curvature term in the action. We find two phases. The weak coupling 
phase is dominated by branched polymer structures with 7 = 1/2. The strong 
coupling phase is flat. The transition between the two phases is first order. We 



show that the existence of the flat phase is not in contradiction with [33|. 

In three dimensions one can interpret the hyper-cubic surfaces also as mem- 
branes. The internal geometry in this model is dynamic, therefore the surfaces 
are called fluid membranes. However, they are not self-avoiding. Self-intersection 
is, differently from typical membranes, allowed. The first order phase transition, 
which we reported in this work, is therefore an example for a phase transition of 
fluid membranes at finite coupling. For fluid membranes described by dynamical 
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triangulations with an extrinsic curvature term the existence of a phase-transition 
flat phase is disputed. In |36] numerical evidence in favor for a second order phase 
transition was reported. The authors of |3^ can not exclude that it is a cross- 
over rather than a phase transition. For fluid membranes with an external curvature 
term in the action a renormalization group analysis ||39|, 40, 41, E2[ indicates that 



there should be no stable flat phase at a finite coupling to the external curvature. 

We cannot take a continuum limit at the transition point because the transi- 
tion is of first order. It is possible that the transition might become softer, and 
continuous, if one introduces self-avoidance into the model. This is for example the 
case for a model of Ising-spins with a gonihedric action [ 43 1 . This model is dual to 
a surface model with extrinsic curvature. It does, however, not define the internal 
connectivity of the surface when more than two plaquettes meet at an external 
link. Numerically a first order phase transition was found [H4[, but with a weak 



self-avoidance imposed evidence for a second order phase transition was found |45]. 
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2.5 The Algorithm 



The algorithm used in the simulations in this chapter is a Monte Carlo algorithms 
which was discussed in section ( |1.7[). I n this section we want to describe the trans- 
formations ( 1.54 ) and the weight (ll.5S| ) specific to our new algorithm for the hyper- 
cubic random surface model. 



The transformations are based on a set of local moves depicted in figure 2.14. 



A prescription for all four pairs of moves is the following: choose p connected 
plaquettes (drawn shaded in the figures) on the surface such that they form a part 
of a sphere (a cube) and replace them by the rest of the sphere, the transparent 
q = 6 — p plaquettes in the figures. We denote by Mpg a move where p is the 
number of plaquettes being removed and q the number of plaquettes inserted into 
the surface. Consider for example the move M15 which removes a single plaquette 
from the surface and fills the hole with a "house" of five plaquettes. In the reverse 
move the role of the shaded and transparent plaquettes is interchanged. For example 
deformation M51 replaces a house with a single plaquette. Deformations M^^a and 
M33fe are self dual, i.e. if one deforms the "same" location twice with one of these 
moves the surface is unchanged. 



Move 15 



s A 

Move 24 



Move 33a 



Move33b 



Figure 2.14: An ergo die set of deformations 



The moves Mpg shown in figure 2.14 are known to be ergodic [p5|. The additional 



moves which we discuss now are therefore not necessary for ergodicity. However, 
they turn out to improve the mobility of the algorithm. The move -/Vfo2 shown in 



figure |2.15| cuts open a link and inserts a spike. A spike are two plaquettes glued 
together along three edges. The move Mqq is a non-local deformation in terms of 
the internal geometry but is local in the embedding space. It cuts open two different 
(internal) links, which occupy the same link in the embedding, and it cross-connects 
the plaquettes. This can probably better be understood in the dual representation. 
In the dual representation a plaquette is replaced by a vertex, links are replaced 
by propagators and a vertex is dual to a closed propagator. However not all closed 
propagator loops are dual to a vertex. One can easily convince oneself that only 
those loops which turn to the right (left) at each vertex of the loop are dual to a 
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Figure 2.15: Spike insertion / deletion 



vertex on the hyper-cubic lattice. The move Mqo cuts open two such propagator 
loops and glues pairwise together the free ends of the two loops, provided this is 
allowed geometrically. In figure 2.16| we show an example of such a transformation. 
The diagram on the left is dual to two spikes glued together along their open links. 
Under the geometric constraint that the two spikes occupy the same unit-square 
in the embedding one can cut open one edge, say the edge to the right from the 
common one, and glue the now free edges pairwise together. After doing this one 
ends up with a "bag in a bag" where a "bag" are two surfaces glued along two 
edges with a common corner. Note that the surfaces may not be glued together on 
opposite edges because the latter configuration has a different topology, that of a 
torus. 





Figure 2.16: An example of move Mqo in the dual representation 

A transformation is done in the following steps: 
1. Choose a location A for the deformation. 



2. Choose a move Mpq from the set of deformations 

3. Choose a direction or orientation d. For example in move M15 one can choose 
in which of the 2(d — 2) possible directions the house should be built. 

4. Check if the move violates any geometrical constraints. If not, subject it to 
a Metropolis test. 

5. Apply the local deformation Mpq at A in direction d. 

The Monte Carlo weights TTpg for deforming the surface at a specific place by move 
Mpq can be derived from the detailed balance condition ( 1.57 ) 



TTpg^ — — e S[A] 
A Um rid 



1 



1 1 



qp 



A + p-qUmu'^ 



re'S[A']. 



(2.28) 



The factor p = ^ is the probability of choosing certain plaquette from the A plaque- 
ttes which form the surface before the move. The grand canonical moves discussed 
above change the size of the surface to A + q — p. The probability of choosing the 
correct location for the inverse transformation is , which is different from p. 
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Therefore the probabihties vTpg contain an area dependent pre-factor to balance the 
grand canonical fluctuations in the transformation probabilities. The next term 
1/nm, the number of possible moves, however, drops out because in the second step 
of the transformation the probability to choose a certain transformation is the same 
for a move and its inverse. The next factor takes into account the freedom one has 
in choosing a direction. For example in move M15 one can choose from 2{D — 2) 
directions orthogonal to the chosen plaquette. In the inverse transformation M51 
one has no freedom after choosing the roof of a house to be removed. The last fac- 
tor is the Boltzmann distribution which we want to generate. Note that the action 
contains contributions from the external curvature. The weights coming from the 
field can be used directly because the curvature field has no degrees of freedom of 
its own. In the following table we list the transition probabilities for the different 
moves. The numbers Ua^ni, give the number of allowed deformations of the type 
Moo at the specific site of the lattice. 



M15/M51 


7ri5 

TTSl 


= 243(d-2)exp(A5) 


M24/M42 


•"'24 


= ^exp(A5) 


M33a/M336 




= exp(A5) 


M02/M20 


•"'02 
•T20 


= 8^(d-l)exp(A5) 


MOO 


TTOO 


= na/rib 



2.5.1 Verification 

To check the validity of the algorithm and its implementation we simulated very 
small surfaces E and compare the distribution of surfaces with a theoretical dis- 
tribution which we have calculated explicitly for the smaller surfaces. To avoid 
too many technicalities we present here in detail only the verification for the model 
without the extrinsic curvature term in the action. In table ( |2.2D we list the smallest 
diagrams. The number N is the number of possible embeddings in the embedding 
space compatible with the internal geometry of the surface. The symmetry factor 
a are also given. Below we will show that the weight of a graph is M ja. We list 
the weight in the columns d = 2, d = 3 for two and three dimensions. In the last 
column the weight measured numerically with our algorithm in d = 3 dimensions 
is given. They are in perfect agreement with the theoretical values. 

In writing the partition function ( |2.2D we were a bit sloppy because we have 
silently dropped the symmetry factor aiE). The justification for doing so is the 
fact that the typical large surface considered so far is not symmetric and hence for 
almost all large surfaces a = 1. The small surfaces we consider now however are 
quite symmetric with a > 1 and hence the factor becomes important. The precise 
partition function is 

Diagrams 

In the rightmost expression we re-express the partition function as a sum over 
diagrams. The number N{E) counts the number of possible ways of embedding a 
surface dual to the diagram E in the Z)-dimensional space. This number is most 
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Diagram 


a{E) 




d=2 


d=3 


N 


a 




8 


Ad{d - 1) 


1 


3 


3.01(2) 


u 





A 


Otl\Uj — ±1 


A 


94 




c 


1 t 

T T 


8 


4d{d - 1) 


1 


3 


3.00(2) 


d 




4 


16d(d - 1)3 


8 


96 


95.97(8) 


e 


o o 


2 


16d{d - 1)3 


16 


192 


192.000 


f 


o 


1 


8d(ci- 1)2 


16 


96 


96.00(9) 


g 




6 


16d{d - 1)3 


5 1/3 


64 


63.99(6) 


h 




12 


4d(d - 1) 


2/3 


2 


2.09(11) 


i 




4 


Ad{d - 1) 


2 


6 


5.97(8) 


j 


(Y) 


24 


8d(cZ- l)(d- 2) 





2 


1.98(10) 


k 




6 


8d(d-l)(d-2) 





8 


8.03(6) 



Table 2.2: Diagrams of size 2,4,6, the symmetry factor a, the number of embeddings 
J\f, the theoretical weight in d = 2 and d = 3 dimensions and the distribution A'" 
obtained numerically in d = 3. The numerical distribution is normalized such that 
diagram e has weight 192. 07 




a 


b 


b 


b 


f 


2 


1 


b 


c 


a 


a 


a 


1 


2 


c 


d 


d 


d 


b 


3 


1 


d 


e 


c 


c 


c 


1 


3 


e 


f 


f 


f 


d 


4 


1 


f 


a 


e 


e 


e 


1 


4 



abcdef 


1234 


++++ 


cdefab 


1342 


++++ 


efabcd 


1423 


++++ 


dcbafe 


1324 


-+++ 


bafedc 


1243 


-+++ 


fedcba 


1432 


-+++ 



Table 2.3: The adjacency matrix (top right) for diagram g using the choice of labels 
given in the top left diagram. The lower table lists the six re-labelings which leave 
the diagram unmodified. 



easily counted in the geometrical representation. In the following we describe how 
to calculate a. 

To depict the diagrams used to verify the algorithm we choose the dual repre- 



sentation in table The hyper-cubic surface model is dual to a M^-theory |4£]. 
The M are hermitian matrices. The graphs are therefore fat graphs discussed in 
context with the matrix model in section ^]^. For clarity, however, we show only 
the thin graphs, known from scalar (f)^ theories, although they do not contain a 
complete description of the geometry. For example a (j)'^ graph does not allow to 
define the geometrical property of two propagators being 'opposite' to each other 
because without the indices carried on a fat graph the propagators interacting at a 
vertex can not always be distinguished. In all graphs shown in this section we use 
the convention to draw propagators which are dual to opposite links on the surface 
as opposite links emerging from a vertex. By propagator loops we denote closed 
propagator loops which go along opposite propagators only, i.e. which do not form 
a right angle at any vertex. Geometrically this is a walk on the surface with one 
internal coordinate kept constant. 

The symmetry factor cr for the diagrams is a discrete remnant of the reparametriza- 
tion group. When drawing a diagram one has a choice of how to label vertices and 
propagator loops. We consider oriented surfaces, therefore one has also a freedom 
in the choice of the orientation of the loops in the dual diagram. Note that all this 
involves only internal properties of the diagram, the embedding is unimportant for 
this. The symmetry factor a counts the number of relabelings which lead to the 
same diagram. For not too large diagrams one can directly probe all possible 

N\ l\ 2' (2.30) 

relabelings for a diagram with N vertices and / propagator loops and count the 



38 



number of invariant relabelings. As an example consider diagram g. In table 
2.3 the adjacency matrix for this diagram with the labels depicted in the left- 
hand figure is written down. One can easily convince oneself that the relabelings 
listed in table 2.3 leave the diagram unmodified. The first transformation is the 
identity transformation, number two has to be interpreted as follows: Rename 
vertex a — > c, 6 ^ d, • • •, loops 1 ^ 1, 2 — > 3, • • • and maintain the orientation for 
loops 1,2,3,4(+ + ++). We found six invariant transformations, the symmetry 
factor for this diagram therefore is o" = 6. 
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Chapter 3 



Four dimensional simplicial 
quantum gravity 

3.1 Introduction 

Quantum gravity in four dimensions, as defined by the path integral ( p,.3[)| ) and the 



action ( f .31 ), is very difficult to deal with. The standard tool perturbation theory 
does not work because the model is not renormalizable. This can be seen from the 
fact that the coupling constant G has a negative mass dimension. Higher order 
loop corrections will therefore always generate more counter terms. 

One possible interpretation of the breakdown of perturbation theory is that gen- 
eral relativity is not applicable at the Planck scale and is just a low energy effective 
theory of some unknown underlying theory. QCD, which is a renormalizable theory, 
provides us with an example of this: the effective low energy theory for the pion is 
non renormalizable. One possible candidate for the underlying theory of quantum 
gravity is string theory. 

On the other hand, a somewhat more conservative interpretation is that the 
divergences in the perturbative expansion do not necessarily mean that general 
relativity does not work at the Planck scale. It may be that they appear just 
because of the naive way we do perturbation theory, we perturb around the wrong 
ground state. The theory may instead be well defined at some non-perturbative 
UV-iiJced. Similarly to the nonlinear a model, which is asymptotically UV free E^], 



the concept of asymptotic safety of four dimensional gravity was introduced |4S]. 
It assumes the existence of a safe sub-space of the coupling space. For couplings on 
the safe surface the renormalization flow leads to a save non-perturbative UV-Hxed 
point. This conjecture is supported by the dimensional 2 + e expansion starting 



from two dimensional quantum gravity, where an UV-&x.ed point was found |4£] for 
e = 2. The model considered in this chapter, four dimensional quantum gravity, 
may provide us with a candidate for the fixed point. 

The structure of four dimensional manifolds is far less understood than in two 
dimensions. For example, there does not exist a topological classification of four- 
manifolds. In this situation the definition of the measure yj^jj jf f^ is very difficult. 



One major conclusion of dynamical triangulations (see section 1.5) is the fact that 
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the sum over discrete triangulations, which replaces ( |1.33| ), defines the measure for 
the two dimensional model in a sensible way. This is a justification for assuming 
that the corresponding sum in the generalized four dimensional model, known as 
four dimensional simplicial quantum gravity, defines a sensible measure as well. 
The assumption implies that the reparametrization degrees of freedom are removed 
in the discrete model. In fact, diffeomorphism invariance is completely lost because 
one can not deform smoothly the discrete structure, which we describe in the next 
section. However, the basic symmetry group of general relativity is reparametriza- 
tion invariance. One therefore has to check whether this invariance is restored in 
the continuum limit of the discrete model, if it exists. 

Simplicial quantum gravity in four dimensions has been extensively studied by 
various groups. The main results are briefly summarized at the end of section 



3.2. In the beginning, the results looked quite promising. In the recent years 
however it became clearer that the model as it stands has some problems. In 
numerical simulations it was found that the model has two phases. The dominating 
geometries in these phases are, however, not extended four dimensional geometries 
as one would expect. The weak coupling phase shows the behavior of branched 
polymers, essentially one dimensional highly branching thin tubes. The strong 
coupling phase is crumpled with a large connectivity. It has an infinite internal 
Hausdorff dimension. The transition between the two phases is first order [ |50| , [5l| ]. 

Based on calculations which assume that mainly the conformal deformations 
are important and ignore transverse excitations of the metric, it was suggested in 
[^2| that the collapse to branched polymers can be prevented by coupling matter- 



fields to the model. This is opposite to what happens in two dimensions where 
additional matter drives the system into a branched polymer phaseQ. In section 3^ 
we investigate the phase diagram of the model with multiple vector fields and show 
that indeed the branched polymer phase disappears. 

Another open question in four dimensional quantum gravity is what is the im- 
portance of topological excitations. It is not a priori clear, whether in the path- 
integral formulation of quantum gravity the sum over metrics should also run over 
topologies. It seems a bit artificial to fix the topology. Let us for the moment con- 
sider spherical topology. One may ask then why space-time can branch but is not 
allowed to join. Branches do not change the topology but if the branches re-join 
the manifold has a hole, i.e. it is not a sphere anymore. In two dimensions it can 
be shown that the naive sum over topologies diverges. In the double scaling limit, 
however, it is possible to define the model with the inclusion of all topological exci- 
tations. Using this two dimensional case as a guidance we have checked numerically 
whether in four dimensions a similar scaling is possible. 



3.2 Simplicial quantum gravity 



The dynamical triangulations model described in section 1.5 is an example how to 
discretize the space of metrics. In this two dimensional model the integral over 
equivalence classes of metrics is replaced by a sum over triangulations. Simplicial 



^See also section 2.1 
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quantum gravity generalizes this discretization to three |^ and four di- 
mensions. In what follows we use D = 3, 4. The integral over equivalence classes 
of D-dimensional metrics is replaced by a sum over D-dimensional simplicial man- 
ifolds. A D-dimensional simplicial manifold is a simplicial complex K with the 
additional constraint that the neighborhood of each vertex p ^ K \s homeomorphic 
to a D-dimensional ball. A simplicial complex is built from d-simplices, where a 
(i-simplex contains d + \ vertices. A 0-simplex is a point, a 1-simplex a link, a 
2-simplex a triangle and 3-simplex a tetrahedron and so on. 

In one and two dimensions a simplicial manifold can be built by gluing together 
d-simplices along their d—\ faces. Each complex constructed in this way satisfies the 
manifold condition. In higher dimensions the resulting complex may however violate 
the manifold condition, i.e. the neighborhood of a vertex is not homeomorphic to 
the D-ball. One can derive relations for the numbers Ni of z-simplices on the 
manifold which are known as Dehn-Sommerville relations [^ ]: 



k=i 



These equations are not independent, they can, however, be used to eliminate the 
A^2i+i in even and N2i in odd dimensions. One can use the discrete expression for 
the Gauss Bonnet theorem (1.32), the Euler relation 



D 



Y^i-iy+^N, = Xd, (3.2) 

i=0 



to eliminate one more Ni, the character Xd is a topological invariant. This means 
that in D = 3,4 only two of the Ni are independent. 

So far we have defined a geometry but not a metric. A Riemannian metric is 
obtained by demanding that: 

• The metric is flat inside a D-simplex 

• The metric is continuous at the D — 1 dimensional faces of the d-simplex. 

• The {D — 1) dimensional faces of a D-simplex S are flat linear subspaces of 
S. 

With this definition the metric inside a L'-simplex S is determined by the length 
of the 1-simplices Si G S. As we will motivate below simplicial quantum gravity 
considers only manifolds with fixed edge lengths l{si) = a. 

To define a path integral of the form ( |1.30| ) for simplicial manifolds we sum over 
all metrics, which can be represented by such manifolds. In the Regge approach 
(see for example ) one considers the metric space obtained by varying the edges 
length ai on a simplicial manifold with fixed connectivity. It is however quite 
difficult to respect the geometrical constraints imposed on the Oj's in the integration. 
Furthermore, a metric is not uniquely parameterized by the Oj's. Therefore one has 
to deal with an unknown weight. Consider for example the triangulation T of a 
flat surface with equal edges length. Moving a vertex in the plane of the surface 
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obviously does not change the metric, although this degree of freedom is integrated 
over in Regge gravity. On the other hand metrics which are very different from the 
manifold T, are represented by singular simplices which contain extremely short or 
long edges length. It is therefore not obvious whether integration over the Oj really 
covers the whole space of metrics. 

Simplicial quantum gravity tries to avoid the mentioned problems by summing 
over dynamical simplicial manifolds T. In other words the dynamics of the model 
now lies in the connectivity of the simplicial manifold. The next step is to pick 
up a parameterization aj = a. This introduces an UV cut-off into the model and 
removes reparametrization degrees of freedom. As in two dimensions one assumes 
that 

V[g]^Y.- (3.3) 

defines the measure in a sensible way. The symmetry factor cr, a discrete remnant of 
the reparametrization group, is important mainly for small manifolds because the 
typical large manifold is not symmetric. In the following we will drop this factor. 

Before one can define gravity using the piecewise linear metric defined above one 
has to discretize the Einstein Hilbert action ( |l.29D . The cosmological term, which 
is the invariant Z)-volume of the manifold is rewritten as the sum of the volumes of 
the D-simplices the manifold is built from : 

d'^eV^^ yD = NDVD. (3.4) 

i^-simplices 

To define the curvature on the discrete lattice one uses the parallel transport of a 
vector on a closed loop on the manifold. In a curved continuous space, the vector 
Va changes by 



^Va = -K,cydfx'dx' (3.5) 

On the discrete manifold one considers a closed loop a dual to a D — 2-subsimplex s 
on the manifold. The loop a connects the centers of the D-simplices Si which have 
s as a common sub-simplex. By definition the space inside the Si is fiat, AT4 = 
on these parts. However when a passes through a, [D — 1) face the orientation of 
Va changes by the deficit angle 

Qd = arccos (3-6) 
in the plane orthogonal to s. The rotation of Sa 

on the closed loop a dual to s is 

therefore described by a single angle 

Rs = 2t:- o{s)Qd, (3.7) 

which we assign to s. The order o(s) is the number of D-simplices which contain 
s as a common sub-simplex. The integrated curvature is a sum over all {D — 2) 
simplices s contained in the manifold : 

^ j dPi^MR ^ 2^ - o{s)eD oc cdNd^2 - ^^^^^^ Nd (3.8) 
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In deriving this we used the fact that a D-simplex has (^jj _ 2) sub-simphces 

and introduced the average order 

27r , , 

of the s in flat space. For example in two dimensions one has C2 = 6, the order of 
all vertices s in the triangulation of a flat surface is o{s) = 6. In higher dimensions 
there does not exist a triangulation of flat space, which is reflected by the fact that 
Cd is non-integer for D > 2. Nonetheless one can ask for the average order of a 
vertex in a hypothetical flat triangulation. 



One gets the discrete Einstein Hilbert action by collecting (3.4) and (|3T 



S = kdNd - kd^2Nd-2, (3.10) 

where the coupling K£)_2 is proportional to the inverse Newton constant G and 
is a linear combination of 1/G and the cosmological constant A. One can consider 



other terms in the action as well. For example in |55] the model with a curvature 
squared term was analyzed, this term, however, does not induce a different critical 
behavior. 

The grand canonical partition function is defined as 

Z{kd-2, Hd) = ^ exp(-Kz)iVD + KD-2ND-2) = ^ eX.p{-KDND)Z{KD-2,ND) 
r Nn 

(3.11) 

where the sum runs over D-dimensional simplicial manifolds T with fixed edges 
length a and fixed topology. The canonical partition function is defined as 

Z{kd-2,Nd) = J2 exp(Kfl_2iVD-2), (3.12) 

where the sum runs over four dimensional simplicial manifolds T/v^ with Nd four- 
simplices and fixed topology. 

We want to close this section with a summary of the results obtained for the 
four dimensional model. The most important question, namely whether or not 



the partition function (3.11) is well defined has been disputed. The exponential 



exp(— K4A'^4) can prevent ( p.llD from diverging only if the number of simplicial 
manifolds M{N4) is exponential bounded. The numerical evidence is in favor of the 



existence of this bound [59, 60, 61|. Recently also an analytic argument |£2[ was 
given for this. The bound defines a critical line ^4(^2) drawn as the solid line in 
the phase diagram figure 3^. For K4 < ^4(^2) eq. ( p. 11 ) diverges and the model 



is undefined. The continuum limit is taken by approaching the critical line from 
above, K4 — > . 

Numerically one finds two phases |5^, |6^. The strong coupling phase is dom- 
inated by highly connected manifolds with very large, possible infinite, internal 
Hausdorff-dimension. The dominating geometries in the weak coupling phase have 
the structure of branched polymers. The internal Hausorff dimension is dh = 2, 
the entropy exponent 7 = ^ in this phase. For the model coupled to matter fields 
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K^-aA +b/GB 



Branched 
Polymer 
Phase 




Crumpled 
Phase X 



Transition Point 
K^- 1.29 



~ 1 / Gb 

Figure 3.1: The phase diagram of four dimensional simphcial quantum gravity 
without matter. 



investigated so far |64, 65, 66] only a renormalization of the coupling-constants was 
observed. 



The strong coupling phase is dominated by two singular vertices |67, 68]. By 
singular we mean that the local volume assigned to these vertices is a finite fraction 
of the total volume in the thermodynamical limit; i.e. the local volume grows 
linear with A'4. The appearance of these singular vertices can be understood in a 
constrained mean-field theory ]7C]. The two singular vertices are joined together by 

2 /3 

a sub-singular link I; its local volume grows like A'^ . The two phases are separated 
by a first order phase transition j^], J^] at K2 ~ 1.29. The singular structure which 
dominates the crumpled phase dissolves at this transition. On small lattices the 
transition actually appears in two steps JTI, 72, ^]. Coming from the crumpled 
phase the sub-singular link disappears first. In the next step the singular vertices 
dissolve. For large volumes these sub-transitions seem to merge [^] . 
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3.3 Topology 



3.3.1 Introduction 

One of the open questions in quantum gravity in dimensions larger than two con- 
cerns the role of topological excitations. It seems a bit artificial to fix the topology 
by hand, which is what is usually done in the simulations. The situation in four 
dimensions is complicated by the fact that as of date no topological classification 
of four manifolds is known. Allowing an arbitrary topology will result in a very 
complicated structure at the Planck scale because curvature fluctuations become 
larger at smaller scales [74|. In [75] it was estimated that the dominant contribution 



to the path integral would come from space-times where the Euler character x is 
of the order of the volume of the space-time in Planck units. Such a configuration 
with many holes is commonly called a space-time foam. 

In 1 76 1 a tensor model similar to the matrix model (section ^]^) was considered, 



which generates pseudo-manifolds with fluctuating topologies. It was shown that 
for large couplings k the pseudo-manifold decomposes in very small connected com- 
ponents while in the weak coupling phase there is only one connected component. 
However, it is not obvious how these pseudo-manifolds, which have now unique 
local internal dimension, are related to four dimensional space-time. 

In this work check numerically, if topology fluctuations are dynamically sup- 
pressed. In a theory containing topology fluctuations only those topologies con- 
tribute to the sum in the thermodynamic limit, which maximize the extensive part 
of the free energy. Other contributions are exponentially suppressed. Two dimen- 
sional quantum gravity provides us with an example of to define a theory which 
contains topological excitations. In the double scaling limit which we discussed 



in section |1.5| , all topological excitations are present. The physical reason is that 
the leading contribution to the free energy is independent of the topology. If one 
believes that all topological excitations are present also in four dimensions the bulk 
volume contribution to the free energy should also be independent of the topology. 
We investigate numerically three different topologies and show that in these cases, 
up to the leading order, the free energy does not depend on the topology. We ob- 
serve a topology dependence in the next to leading order. We study those volume 
corrections and analyze their kinematic sources. 

3.3.2 Model and methods 

Given the canonical partition function ( 3.12| ) , we define the intensive free energy as 



/(/t2)= lim -^logZ(AC2,iV4). (3.13) 
For large A'4 the free energy is assumed to have the form: 

F(/t2,A4) = A4/(/t2)+5(A^2,A4), (3.14) 

where the function (5 is a finite size correction, i.e. for any K2 

hm i%^ = o. (3.15) 
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The derivatives of the free energy are convenient to recover some basic properties 
of the theory. The action density 



is normahzed such that it becomes an intensive quantity in the thermodynamic 
hmit. We call the derivative with respect to 

dF 

'^4(«2,iV4) = ^. (3.17) 

the pseudo-critical value of the parameter K4. This definition differs slightly from 
the one proposed in |5£]. This quantity is a measure of finite size dependence of the 



free energy, i.e. of how fast it approaches the thermodynamic limit K4 — > K4{k2) 
for ^00. In this limit the value ^4(^2) defines the critical line of the model 



corresponding to the radius of convergence of the series ( 3.11 ) 



The pseudo-critical value ^4(^2, A^4) is related to the action density by 

Oka dr , . 

where the second term goes to zero when A'4 goes to infinity. Therefore the critical 
parameter K4 becomes an integral of r in the thermodynamic limit. Thus, if r is 
independent of the topology so is K4 unless the integration constant depends on the 
topology. To fix the integration constant it is sufficient to measure K4 for only one 
particular value of K2 ■ 

The action density ( 3.161 ), the normalized average number of triangles on the 



manifold, is a standard observable which is easily accessible in canonical simulations. 
For the pseudo-critical value K4 ( 3.17| ) more unconventional methods are required. 



The baby universe method, which we used in the context with hyper-cubic random 
surfaces to extract the next to leading correction, unfortunately works only in the 
elongated phase. Following |59|, we adopt here a more general method based on 



grand-canonical simulations, which works equally well in all phases. We use multi- 
canonical simulations with a potential U{N4). From the volume fluctuation one 
can extract the volume dependence of K4. We expect the quantity dF/dN4, which 
we want to measure, to smoothly approach the function /(K2). The freedom one 
has in choosing the potential U can be used to minimize the statistical error of the 
observable. A Gaussian potential controlling fluctuations around a fixed volume V4 
is well suited for this problem ||5^ : 



U = K4N4 + ^{N4 - V4f (3.19) 



but other potentials can be used as well |55]. With this potential the multi-canonical 
partition function reads : 

Z = gK2iV2-(7(A^4) ^ ^ gF(K2,Af4)-K4Af4-^(iV4-V4)2 {3.20) 

T N4 
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Figure 3.2: A four tetrahedron built into a four dimensional hyper-cube. 



The distribution of the volume A'4 becomes Gaussian around A'4 = V4 if one expands 
F(k2, A^4) around this point 



P{x = iV4 - Vi) ~ exp(--(x - xof + ...) 
The parameters T and xq are related to F and U by: 

r = j-d^F/dVi, xo = {dF/dVi- Ki)/r. 



(3.21) 
(3.22) 



Both quantities can be measured in the simulations, by tuning K4 close to the 
derivative dF/dV^ and choosing 7 much larger than the second derivative d'^F/dV^. 
The parameter F can be extracted from the width of the distribution and xq from 
the shift of the maximum from A'4 = V4. This gives an estimate for the pseudo- 
critical coupling at V4 : 

K4 = dF/dV4 = Txo + Ki. (3.23) 

In this formula, K4 is the coupling used in the simulation from which xq and F were 
extracted. The estimator can be improved by minimizing xq. This can be done 
recursively by replacing K4 by K4 in subsequent simulations. 

The value for the multi-canonical coupling 7 should be chosen large enough so 
that the volume fluctuations are small compared to the scale of change of F{N4, K2). 
On the other hand, the algorithm used in the simulation is grand canonical and 
requires volume fluctuations for ergodicity. Therefore 7 should be small enough to 
allow sufficient A'4 fluctuations. We checked that for 7 = 0.0001, F computed from 
the width of the distribution F = a'^^N^)^ i.e. AA'^4 = 100, was equal to 7 
within the error bars. From (|3.22|) it follows that 7 > d'^F/d'^VA. Therefore the 
free energy changes very slowly across the range of the A'^4 distribution width as 
needed for the approximation ( ^.21 ) to be valid. 

Because the local moves fl% used in the Monte Carlo simulation preserve the 
topology t of the manifold, t is given by the topology of the starting configuration. 
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In our simulation we simulated a spherical topology and two different tori, 
namely x and x x x S\ 

The starting configuration for the sphere is, as usual, the 4d-boundary of a 5- 
simplex. The x x x manifold can be built out of the regular 3^ square 
torus by dividing each elementary 4d-cube into 4-simplices in the following way. 
For each cube we mark two points, pi = (0, 0, 0, 0) and = (1, 1, 1, 1), lying on the 
opposite ends of the main diagonal and connect them by one of the shortest paths 
going along the edges of the cube. The shortest path goes through 4 edges, each 
in a different direction, and through three points, say pi,P2,P3- There are 24 such 
paths. Each set of points P0iPiiP2,P3,P4: forms a 4-simplex. In figure 3^ we show 
as an example one four-simplex (red) built into a four dimensional cube obtained 
with this procedure. The dashed lines are additional diagonals which have to be 
inserted to complete the four-tetrahedron. The manifold condition requires that 
no vertex be connected to itself. Therefore one has to glue together 4 cubes in all 
four directions. The starting configuration therefore contains 3^ * 24 = 1944 four 
simplices. 

The manifold x can be produced from a spherical manifold by taking 
away two separate 4-simplices and gluing together the boundaries bi , 62 created by 
removing the two 4-simplices. However, two vertices in a 4d simplicial manifold 
may have at most one common link. To avoid the creation of a double connection 
in the gluing process, the distance between the vertices on bi to those on 62 has to 
be at least three links. One can ensure this, without going through a tedious check, 
by gluing cyclically together three copies of the original (spherical) manifold, i.e. 
bi/b2,b'i/b2,b'(/b2. The (double-) prime is used to distinguish the different copies. 
The boundaries 61, 62 are chosen such that they have no common vertex. 



3.3.3 Numerical results 



We have performed simulations with N4 = 4000, 8000, 16000, 32000, 64000 for the 
three topologies discussed above. The explicit form of the sub-leading correction 
(5(k2,A^4) depends on the phase where the measurement is taken. The elongated 



phase of simplicial gravity is dominated by branched polymers [55, 71, 78]. This 
means that in the large volume limit one can use the ansatz 



F(K2,iV4) 



+ (7 - 3) log7V4 + M^2) + O {^^^ 



(3.24) 



for the free energy. The correction coefficient 7 is assumed to depend only on 
the genus g of the manifold. The logarithmic corrections disappear for the action 
density r eq. ( |3.16| ) 



1 dF{K2,Ni) 
Ni 8^2 



/o(^2) + 



N4, 



(3.25) 



Therefore one should take the next order corrections 1/A^4 into account. We will 
show that they appear for purely kinematic reasons. Consider the limit of large pos- 
itive K2 in which only manifolds maximizing N2 contribute to the sum (PI ). Such 
triangulations can be obtained from barycentric subdivisions (3.54) of 4-simplices 
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applied successively to a minimal starting configuration, possibly mixed with micro- 



canonical transformations (3.74), which do not change N2 and A^4. By the minimal 
configuration we mean the minimal volume triangulation which maximizes For 
the barycentric subdivisions one gets the relation = 5/2A^4 + c^, where the 
constant 

cO = ivO _ 5/2ivO (3.26) 

characterizes the initial minimal configuration. The number of triangles A^2 is 
related to the action density r = A^2/-^4 = 5/2 + /N/^. This means that the 
constant leads to I/N4 corrections of the action density. The contributions to 
the sum ( p. 11 ) from triangulations built from non-minimal ones (i.e. smaller 



are suppressed exponentially by exp— ^2(0'^ — c), where c characterizes the non 
minimal starting configuration with 0(^2) < c^. 

The minimal configuration for the sphere, the surface of a 5-simplex, is known 
and hence is known. For the two toroidal geometries we find the minimal con- 
figuration numerically. We use our Monte Carlo code in a cooling procedure, in 
which we increase K4, to decrease and increase K2, to maximize A^2- In practice 
we have to increase K2 slowly compared to K4, because increasing K2 also increases 
pseudo-critical coupling K4{k2). For the topologies we studied, we found the fol- 
lowing minimal configurations : 



(3.27) 







n2 


= 6, 


Nl 


= 20 




= 5 




X 


K 


= 110, 


Nl 


= 44 




= 




xS^ 


n2 


= 704, 




= 1472 




= 2 



For the manifold 5"^ x the finite volume correction l/A'4 is not present at all. 
It is also very difficult to detect them for the sphere because c'' = 5 is very small 
compared to the volume A'4 used in the simulation. For the torus x x x 
the corrections are, however, two orders of magnitude larger and measurable in the 
volume range used in the simulations. This estimate of the 1/A^4 effect is exact only 
for infinite positive K2 but one expects it to work, although with a slowly varying 
coefficient c(k2), in the entire elongated phase. 

In figure |3.3| we show the action density r measured in the elongated phase for 
K2 = 2.0 against l/A'4. In the same figure we display the curve Too -|- (P /N^ with 
= 288 and = 2.482 which fits the data points very well. Note that that r^o 
does not depend on the topology. 

For the pseudo-critical coupling K4 we observe, with the statistics available, no 
volume and/or topology dependence either. We find kf = 5.659(4) for K2 = 2.0 in 
the infinite volume limit. This is compatible with the ansatz 

^4(^2, A^4) = /o(^2) + (3.28) 

because 7 is known to be of order 0(1). With the method used, the correction 
could only be separated from the statistical noise with an disproportionate amount 
of computer time. 
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IOOO/N4 

Figure 3.3: The volume dependence of the action density r for three different 
topologies in the elongated phase at K2 = 2.0 

In the crumpled phase we use a power-law ansatz 

F{k2, N^) = N^U^2) + h{^2)Ni (3.29) 
Taking the derivative with respect to one gets: 

ki{K2,Ni) = h{K2) + 5h{K2)Nl-\ (3.30) 
For the action density ( ^.16|) one finds 

r = at^2) + f[{^^2)Nt^ (3.31) 
We checked this ansatz by fitting our numerical data for r{N4) to eq. ( ^.311) : 



Topology 


6 




log(/() 


xVdof 




0.5 (2) 
0.6 (2) 
0.5 (2) 


2.039 (+0.010 / -0.013) 
2.028 (+0.008 / -0.016) 
2.038 (+0.010 / -0.021) 


1.9 (1.3) 
1.0 (0.9) 
1.8 (1.3) 


0.78 
0.87 
0.31 



and K4(iV4) to (|3!30| ): 



Topology 


6 


1,00 f 

1^4: — JO 


log(5/i) 


xVdof 


X 


0.6 (2) 
0.6 (3) 
0.6 (2) 


1.20 (2) 
1.20 (2) 
1.20 (2) 


1.4 (1.3) 
1.4 (1.5) 
1.4 (1.3) 


0.22 
0.10 
0.18 



In the tables we give the logarithms of /( and 6fi since these have approximately 
symmetric errors. One can see that the values do not depend, within errors, on 
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1.06 



0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 



Figure 3.4: The volume dependence of the critical «;4 for three different topologies 
in the crumpled phase at K2 = 0.0 



topology. In figure 3.4 we show the numerical data for K4 for the three topologies 
and 5 = 0.5. 

The source of the volume correction in this phase can be understood by looking 
at the typical configurations which dominate the ensemble for large negative K2- 
These configurations minimize the free energy and therefore have for fixed volume 
the minimal number of vertices, and thus also the minimal number of triangles. In 



three dimensions such configurations were constructed in [57|. The generalization 
to the four dimensional case is obvious. One starts with a 2d triangulation of the 
sphere with t triangles. At each triangle one builds a Ad-pancake neighborhood 
from q 4-simplices lying around the triangle in such a way that the links opposite 
to this triangle form a circle. An opposite link is defined as a link which does 
not share a common vertex with the triangle. The next step is to put neighboring 
pancakes together by identifying these circles with the 3d faces of the neighboring 
pancakes. Each pancake has three such faces which lie between the circle and an 
edge of the basic triangle. It also has three neighboring pancakes. After this step 
one gets an sphere with t * q 4-simplices and 2 + + g vertices. The highest 
connectivity ~ {Nq — 2)^/2 is reached when t = 2q. The number N2 of triangles 
can be constructed from this as a function of A'4 making use of the Dehn-Somerville 
relations ( |3.1| ). Inserting the result into ( p.l6| ) one finds r oc 1/ / \/\N/i) i.e. 5 = 0.5 
for these configurations. 

In 1^^, evidence was given that for spherical topology the phase transition is 
of first order. We could observe flip-flops in the action density for the two tori as 
well. This seems to show that the transition in these cases is of first order as well. 
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3.3.4 Discussion 



We have investigated the behavior of the curvature and the entropy density for 
three different topologies in four dimensional simplicial quantum gravity. We con- 
centrated on two K2 values, one deep in the crumpled phase and another value in 
the elongated phase. In both cases the value of the entropy density ^4(^2) and the 
curvature were independent on the topology in the large volume limit. We have 
shown that is an integral of the curvature. The observation that k% is indepen- 
dent from the topology in the two phases of the model therefore suggests that the 
curvature and the entropy density do not depend on the topology for all K2 in the 
thermodynamical limit. This gives support to the conjecture that it is possible to 
define, similar to the double scaling limit in two dimensions, a scaling with coherent 
contributions from all topological excitations. We can, however, not fully exclude 
that in the transition region the curvature and hence K4 depend on the topology 
but this, by chance, cancels out after integrating over this region; this seems very 
unlikely, though. 

Our result furthermore supports the conjecture that the entropic bound exists. 



This question was recently discussed in the literature 8C, 81 1. The value of 
the entropy density K4 can be compared with the estimates given in [62|, based 
on a summation over curvature assignments. This estimate, does, however, not 
directly give results in agreement with our numerical data for K4. We conclude 
that the mapping from the entropy of curvature assignments to the entropy of four- 
manifolds, which is a basis for the calculations in [62|, is to simple to capture fully 
the behavior in the AD approach. 

We have analyzed the finite size effects observed for the model. We gave kine- 
matic arguments which explain the leading correction very well. 
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3.4 Simplicial quantum gravity coupled to gauge mat- 
ter fields 



3.4.1 Introduction 

As mentioned in the introduction, the dominating geometries in the two phases 
found in four dimensional simplicial quantum gravity do not describe extended 
four dimensional structures. This, however, is required for the model to describe 
the large distance behavior of space-time. The phase transition is of first order. 
Therefore one can not take the continuum limit at this point. In this situation it 
is sensible trying to insert terms into the action which cure these problems. 

In the crumpled phase, for example, one may try to suppress configurations 
with large order o(t), i.e. singular vertices, by adding a term 



alog^o(t) (3.32) 



to the action. We have simulated the model with this term. We found, however, 
almost no back-reaction of the gravity sector on this term. Both phases survive. 
At the transition one finds flip-flops in the time-series of A^o which indicates that 
the transition is still flrst order. 

The branched polymer behavior observed in the elongated phase in fact plagues 
many theories of random geometries in all dimensions. One phase of the hyper- 
cubic random surface model investigated in this thesis is a example of such a theory. 
On the other hand, pure gravity in D = 2 dimensions is in the non-trivial Liouville 
phase with string susceptibility exponent 7 < 0. In this phase it is possible to deflne 
a sensible continuum limit. The branched polymer behavior returns, however, if 
one couples matter with central charge c > 1 to two dimensional quantum gravity. 



In |82, |8^, 84] a heuristic argument was given to explain this behavior: When 



c becomes larger than the C), = 1 'barrier' metric singularities condense. The 



argument is similar to the Berezenskii-Kosterlitz-Thouless |85, I83] argument for 
a transition in the 2-D x — y model. It relies on the construction of singular 
conflgurations ('spikes') and the observation that the action for such conflgurations 
grows logarithmically with the infrared cut-off. From the continuum methods and 
the known dynamics of the conformal factor ( |1.23 ) it is known that the action for 
a spike 

5(c) oc (26-c)log(-) (3.33) 
a 

becomes smaller with larger c. The entropy E oc log(-L/a) of putting spike onto 
the surface on the other hand grows logarithmically with the linear extent L of the 
surface. It was shown that for c > 1 the entropy E becomes larger than the action 
5(c) and the creation of spikes is favored. 



In 1 87] it was suggested that a similar phenomenon might occur in four dimen- 
sions. An effective action for the conformal factor in four dimensions was derived 
in 



89]. As in two dimensions the action for a spike grows logarithmically with 



the linear extent L : 

SeWs] = \Q'q' log (^]. (3.34) 
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The couphng 



= + Y^^^ + 62iVy - 28) + Ql,,, (3.35) 

depends on the numbers A''^ of scalar fields, Nyyp of Weyl fermions and Ny vector 
fields. The term Q'^^av describes the (unknown) contribution of spin-2 gravitons. 
The only negative contribution —28 comes from the conformal factor itself. The 
most interesting point about this result is that, differently from the two dimensional 
case, the effective action becomes larger when one adds matter to the model. In 
other words, putting matter into the model stabilizes the model and suppresses 
spikes. This suggests that in the absence of matter fields quantum gravity has no 
sensible vacuum, but this instability is lifted when the number of matter fields is 
large enough. 

It should be emphasized that in four dimensions the action ( 3.34| ) for spikes 



can not be derived rigorously, because it is not clear how transverse excitations 
of the metric should be treated. Nevertheless the conjecture above makes is very 
interesting to investigate the influence of matter on the geometry in 4d simplicial 
quantum gravity. It may also explain why in earlier investigations ^ no 

non-trivial effects induced by the matter have been observed. In the framework of 
[^2| it was shown that spikes are dominant if < S/q^. It may be that in the 
earlier investigations not enough matter fields were introduced so that was still 
too small. 

3.4.2 The model 



From ( 3.35 ) one can read off that the strongest effect on the geometry is caused by 
the number of vector fields. Therefore we study a model of continuous vector gauge 
fields coupled to 4d simplicial gravity. Specifically we consider non-compact U{1) 
gauge fields A{lab) living on the links lab of the simplicial manifold and A{lab) = 
—A{lba)- The sum in the action 

Sm = Y1 0{tabc) [A{lab) + A{lbc) + A{l,a)f (3.36) 
^abc 

extends over all triangles tabc- The order o{tabc) is equal to the number of four 
simplices which contain t^fcc and is proportional to the local volume assigned to the 
triangle. This factor appears to mimic the factor which appears in a covariant 
action. There is no need to introduce a coupling in front of the sum because the 
action has a Gaussian form. 

In the model with n gauge fields we put n replicas of the matter field on the 
dynamic geometry. The action 



= E ^'^Sm (3.37) 



is the sum of the one field actions (|3.36| ) with ^^^Sm = Sm- Note that there is no 
explicit interaction between the replicas of the model. On a fixed lattice such a 
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model is trivial, the behavior of a single replica remains unchanged. On a dynamic 
lattice, however, the systems are coupled by the back-reaction of the dynamical 
geometry on the matter fields. 

The action for the combined system of gravity and matter contains, besides 

(n) 

Sj^ , two more terms. The contribution of the gravity sector is given by the discrete 
Einstein-Hilbert action (|3.1[l| ) which in four dimensions reads 

Sg = -K2N2 + K4N4. (3.38) 

where A''^ denotes the number of /c-simplices. The third term is present for tech- 
nical reasons. The algorithm in the simulation is inherently grand canonical and 
requires volume fluctuations to satisfy ergodicity. We use the multi-canonical poten- 
tial ( |3.19| ) to simulate a pseudo-canonical ensemble with almost fixed ~ 

Measurements are taken only if A'4 = A'^^'^^ The coupling K4 is tuned so that 
JU = N4. 

The model is defined by the partition function 

T '' i=l \leT / 

(3.39) 

The sum extends over all distinct four dimensional simplicial manifolds T and W{T) 
is the symmetry factor taking care of equivalent re-labelings of vertices. The prime 
at the integral indicates that the zero modes of the gauge field are not integrated. 



3.4.3 The strong coupling expansion 

The first few terms of the strong coupling expansion can be used to estimate the 
exponent 7 |9^. We used our code without matter fields to identify diagrams up 
to A'4 = 30 and calculated the symmetry factor a for the diagrams by counting 
the number of invariant re-labelings 0. Without matter the weight for a diagram is 
W = 1/a. To debug our computer code we identified the diagrams with A'4 < 18 
by inspection and verified that the theoretical and numerical distributions governed 
by the symmetry factor agree. 

For the strong coupling expansion of the model with matter-field we have to cal- 
culate the weight induced by the matter. The Gaussian integration which excludes 
zero-modes over one species of fields leads to a determinant p. As the n replicas 
of the field interact only indirect via the back-reaction of the gravity and have no 
explicit interaction term in the action it is sufficient to calculate the weights pM for 
one field. The n-field weight is 

= pl^/a. (3.40) 



The weights pM were calculated using Maple. In table ( |3.lD we summarize the 
results obtained with this procedure. 

The series was then analyzed with the ratio method |90| to extract the exponent 



7. In the upper plot in figure |3.5| we show the behavior of 7 as a function K2 



For a detailed discussion of this method see section 2.5.1 
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No 


^9 














6 
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1 


1 


1 








10 


7 


1 


3 


0.097638467.. 


1.03423.. 


. X 


10- 


-4 


12 


7 


1 


5 


0.030058406 . . 


1.08632 . . 


. X 


10" 


-6 


14 


8 


1 


15 


0.018550484.. 


2.83715.. 


. X 


10^ 


-8 


16 


8 


2 


255/4 


0.015777808 . . 


9.68553 . . 


. X 


10^ 


-10 


18 


8 


3 


110 


0.005500465 . . 


1.38182.. 


. X 


10- 


-11 




9 


3 


95 


0.004759295 . . 


1.19996.. 


. X 


10- 


-11 


20 


8 


2 


225 


0.002512817.. 


3.15034. . 


. X 


10- 


-13 




9 


7 


693 


0.007291315.. 


8.14758.. 


. X 


10- 


-13 


22 


9 


15 


2460 


0.005728290 . . 


3.27573 . . 


. X 


10- 


-14 




10 


7 


690 


0.001447804 . . 


6.46761 . . 


. X 


10- 


-15 


24 


9 


13 


16365/2 


0.004226212.. 


1.17586.. 


. X 


10- 


-15 




10 


34 


14625/2 


0.003378959 . . 


7.45244 . . 


. X 


10- 


-16 


26 


9 


50 
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0.001946262 . . 


2.34367.. 


. X 


10- 


-17 




10 
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39645 


0.003936950 . . 


3.97116. . 


. X 


10- 


-17 




11 


30 


5481 


0.000491334.. 


4.06700 . . 


. X 


10- 


-18 


28 


9 


89 


291555/7 


0.001058334.. 


6.96159.. 


. X 


10- 


-19 




10 


415 


182820 


0.004119603.. 


2.16619.. 


. X 


10- 


-18 




11 
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77057 


0.001534637.. 


6.33995 . . 


. X 


10- 


-19 


30 


9 


139 


73860 


0.000457581 . . 


1.78973.. 


. X 


10- 


-20 




10 


1276 


672821 


0.003427165 . . 


9.13722.. 


. X 


10- 


-20 




11 


1208 


564000 


0.002507268 . . 


5.17157.. 


. X 


10- 


-20 




12 


143 


46376 


0.000179907.. 


2.84432 . . 


. X 


10- 


-21 



Table 3.1: The number of different graphs Ng, for a fixed volume A^4 and fixed 
number of vertices Nq and the corresponding weights WNyiN^, Nq) for 0, 1 and 3 
vector fields. All weights are normalized with the value at A^4 = 6. 
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Figure 3.5: (a) Variations of 7 with for 0, 1, 3, and 5 vector fields coupled to 
gravity, (h) Variation of 7 with the number if fields for K2 = 10.0. These values 
are obtained using the ratio method to analyze the strong coupling scries including 
terms corresponding to A^4 = 6, 10, 14, 18, 22, 26, and 30. The horizontal line in 
the lower diagram is at 7 = 0.5, the exponent expected for pure gravity. 
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Figure 3.6: The phase diagram for the model with one vector field. 



obtained with this method. For large K2 the susceptibility exponent 7 tends to a 
constant value which becomes negative when the number of fields grows. For pure 
gravity one finds 7 = ^ as expected. For small enough K2 the exponent 7 is large 
negative. In the intermediate region the results are unstable. This presumably 
indicates a phase transition in this range of the coupling K2- In the lower plot we 
show the variation of 7 as a function of the number of vector fields at K2 = 10.0, 
i.e. in the weak coupling phase. We used diagrams of size [6, . . . , 22], [6, . . . , 26] 
and [6, . . . , 30] for the three different curves. One sees that they depend slightly on 
the size of the largest diagrams used in the series expansion, however they seem to 
converge for larger diagrams. This and the numerical results presented in the next 
section justify the assumption that the series expansion presented here describes 
correctly the qualitative behavior of the exponent 7. 



3.4.4 Numerical results 

The first question when analyzing a new model is about the phase diagram. The 
gauge fields investigated here introduce a new parameter into the model, namely 
the number n of fields. We have performed Monte Carlo simulations for the model 
with one and three vector fields. For the survey of the phase structure we use rather 
modest lattice sizes with no more than 16K four-simplices. It is known from the 
simulations of pure gravity that finite size effects are quite strong; for example the 
order of the phase transition for pure gravity reveals itself only for lattices larger 
than 32K. It is therefore clear that one has to be careful when drawing strong 
conclusions from the numerical data presented here. Nevertheless we believe that 
the qualitative picture obtained from the data is reliable. This is supported by 
the observation that the results are in good agreement with the strong coupling 
expansion. 



In figure 3.6 we present the phase diagram found for one vector field. The model 
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diverges if K4 < ^^(^2) and the model is undefined. We observed no back-reaction in 
the strong couphng phase, which is still dominated by singular structures. Likewise, 
the weak coupling phase is still a branched polymer. In the transition region we 
observe two sub-transitions indicated by two peaks in the node susceptibility 

Xo = {{Nl) - {Nof)/N^. (3.41) 

The first sub-transition has a broad peak. With the statistical noise it is difficult 
to locate precisely the point of the transition. Therefore our data does not allow 
to judge if with growing volume the two transitions merge to the same critical 
coupling. They are presumably similar to the two sub-transitions observed for pure 
gravity on small lattices^, which do seem to merge |7^ . 

The picture changes when we include three vector fields. We still find singular 
structures in the strong coupling phase, but the weak coupling phase has no longer 
a branched polymer behavior. The two phases are separated by a maximum of the 



node susceptibility at K2 « 1.9 which grow with the volume as seen in figure 3.7. To 
locate the critical coupling and understand the nature of the transition one needs 
better statistics and simulations for larger systems; this is in progress but it will, 
however, require several more months of CPU time. 

The change in the geometry at roughly the same coupling where xo has its 
maximum is, however, a strong indication that we observe a real transition. In plot 



(b) of figure 3^ we show the difference between the orders of the most singular and 
the third most singular vertex pi — ps- In the strong coupling phase the order of 
the first and second most singular vertices are roughly of the same order. They 
are seperated by a large gap from the order of the other vertices. For larger K2 the 
order of these two vertices merge to the order of other vertices. 

This means that there is no singular vertex in the weak coupling phase, as 



demonstrated by the scaling behaviour of the most singular vertex pQ in figure 3.5 
for K2 = 4.5. The local volume of this vertex grows sublinearly like N^^'^, i.e. a 
vanishing fraction of the total volume in the thermodynamical limit. 

In the same plot (13.81) we show the scaling of the average distance between the 



simplices [91|: 



were n(r) counts the number of simplices at a geodesic distance r from a marked 
simplex. With a fit to equation ( 3.42| ) we get for the Hausdorff dimension = 



3.93(15). However, the data points do not lie on a straight line, presumably due to 
finite size effects. If we exclude the data point for the smallest volume = 2000 
we get dh = 4.68(8). 

We use the baby universe method |9^] to measure the entropy exponent 7 in 
the weak coupling phase. The method was explained for the hyper-cubic random 



surface model in section 2.3.2. Up to some technicalities the method is the same for 



the model discussed here. In table 3.2 we list the results obtained for three vector 



fields with this method for two values of K2 (4.5 and 6.0) and at different volumes. 



See also the end of section p. 2 
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Figure 3.7: (a) The node susceptibility xo, tor three veetor fields coupled to gravity, vs. 
K2- For = 16000 the statistics is not sufficient for a reliable estimate of xo- (b) The 
corresponding change in the difference between the orders of the first and of the third most 
singular vertex, po and P2, normalized with the total volume N^. 





K2 = 4.5 


K2 = 6.0 


2000 


-0.22(2) 




4000 


-0.18(3) 


-0.17(4) 


8000 


-0.23(3) 




16000 


-0.30(6) 


-0.12(6) 



Table 3.2: Measured values of the string susceptibility exponent 7 in the weak 
coupling (large ^2) phase for three vector fields 
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Figure 3.8: The scaling of both the largest vertex order po ^^nd the average distance 
between two simplices in the new weak coupling phase, for n = 3 vector fields 
at K2 = 4.5 in a doubly logarithmic plot. The lines are the best linear fits to the 
data. 

We consistently observe negative values 7 ~ —0.2. This is in reasonable agreement 
with the prediction from the strong coupling -0.38. 

3.4.5 Discussion 

The most interesting result is the discovery that matter fields can prevent the 
collapse of 4d simplicial manifolds to branched polymers. We find for the first time 
a strong back reaction of the matter on the geometry. The results suggest that 
when the number of vector-fields becomes larger than one the branched polymer 
phase disappears, and is replaced by a new phase with negative 7. 

The numerical data does not allow to determine the order of the transition 
unambiguously. However for one vector field we observed indications for a first order 
transition but the latent heat of the transition is smaller compared to the jump 
observed for pure gravity. This may indicate that matter softens the transition. 
Therefore it is not excluded that with the correct mixture of matter-fields the 
transition is continuous. 

The new weak coupling phase observed for three vector fields is very intriguing. 
The internal Hausdorff dimension if close to four, the same as for flat space. The 
negative 7 values obtained numerically and with the strong coupling expansion 
indicate a non-trivial behavior. This is a very promising result because the negative 
7 may, like in two dimensions, indicate that the whole phase is critical and a 
non-trivial continuum limit can be taken. Gravitons, however, will presumably 
be present only if the Newton coupling is tuned to a critical value k^- If the 
transition found in our work is second order this induces infinite curvature-curvature 
correlations, which may be interpreted as gravitons. 
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Our results obtained numerically and with the series expansion support the con- 
jecture in 1 52] mentioned in the introduction. As predicted we find that additional 
matter stabilizes the geometry as opposed to the case in two dimensions where 
matter de-stabilizes the surface. 

After all it should, however, not be forgotten that one still faces some uncer- 
tainties. It is important to check carefully the nature of the transition. It should 
also be verified that the critical coupling in the thermodynamical limit is finite. 
In the weak coupling phase the curvature N0/N4 « 0.25 is close to the kinematic 
bound. The local volume o{p) oc N^^^ of the sub-singular vertices in this phase 
take no finite fraction of the whole volume in the thermodynamical limit, but the 
local volume still diverges. 
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3.5 Algorithm 



The algorithm used to generate the data presented in this chapter uses the frame- 
work of dynamical Monte Carlo algorithms discussed in section (|1.7D . In this section 



we want to describe the transformations ( 1.54 ) and the weight (|l.59D specific to the 



algorithm used to simulate four dimensional simplicial quantum gravity. We also 



discuss the weights induced by the vector field used in section 3.4. 

In order to calculate the transition probabilities, remember that the distribution 
of manifolds T we want to generate is given by 

p{T, {A}t) ~ exp -Sg[T] - Sa[T, {A}t], (3.43) 

where the first term corresponds to the Einstein-Hilbert action. The action Sa of 
the matter field living on the manifold is given by 

n 

Sa[T,{A}t] = ^^o(0(^(^)(/i2) + ^«(/23) (3.44) 

i=l t£T 

The second sum runs over all triangles t of the triangulation T, o{t) is the order 
of the triangle t, and /12, ^23 1^31 are the three oriented edges of tl23. The field A 
is an oriented link variable, which means A{—1) = —A{1). Therefore the action 
is invariant under a gauge transformation at vertex pq: simultaneous translation 
A{1) A{1) + e of all A{1) living on links emerging from po does not change the 
action. The first sum runs over the number of matter fields coupled to gravity. 
For the sake of simplicity, we discuss in the following only the case n = 1. The 
extension to n > 1 is straightforward. 

The geometric moves M(^p^q-j which we describe in detail below may generate 
and destroy links and hence field variables A{i). In principle one could update the 
geometry and the matter at the same time. A procedure would be to assign random 
values to new matter variables and accept or reject the whole transformation with 
a Metropolis weight. However, if one assigns the A{li) with a flat distribution 
the efficiency is very low because then most transformations will be rejected in the 
Metropolis step. One might try to assign to the configuration A a Gibbs distribution 

V{A) = [ P^e-^-^(^s) (3 46) 

J A 

The integral runs over all accessible states of the new matter-fields. The prime 
indicates that we do not integrate over the zero-modes of the action. The prob- 
ability can, however, not directly be used as a heat-bath probability because in 
the grand canonical moves the normalization constant V{A) is different for a move 
and its inverse. Therefore our algorithm updates gravity and matter in two steps. 
In the first step the geometric part of the transformation is accepted or rejected 
independently of the state of the new matter fields, which are integrated out. In 



with the normalization 
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the second step the fields which were created in the geometric update are set with 
a heat-bath probability. 

The detailed balance condition for the probability ^(p^q) to accept the geometric 
move {p, q) at a certain location on the manifold A reads then 



ViA) 



n(p,q)(^)e' 



-Sg(A)-Ss{A) 



^n(,,)(^)e-5«(^)-^^(^). (3.47) 



The term Sg is the action of the matter field without contributions from terms 
involving the new matter field. In the geometrical update Sg changes because the 
move changes the order o{t) of some triangles on the simplicial manifold. The 
factor involving the rii are present to balance the grand canonical fluctuations of 



the manifold. The n,- 



3, 4, 5 count the number of possible locations, where a 



local deformation can be applied. We have discussed this factor in detail in the 
context with expression ( 2.2^ ) for the hyper-cubic random surface model. 

The matter field is not only updated in the geometrical transformations. The 
algorithm switches between updates in the geometrical and the matter sector. In 
the matter update we choose an extensive number of links from the manifold. Each 
link is sequentially subjected to a heat-bath or an over-relaxation step. In the 
heat -bath, link field A on the link I12 is chosen from the distribution : 



P{A) ~ exp ^- o{t){A - Alp + A2pf^ = exp {-a{A - 13/af 



(3.48) 



where the sum runs over all triangles t = ti2p which have the link I12 as an edge and 
the point j? as a vertex. The variables a, (3 in the rightmost expression are defined 
by 

a = 'Yo{t) (3.49) 



P = Y^o{t){Aip-A 



2p>- 

The over-relaxation is a micro-canonical transformation 

A -A + 2f3/a, 



(3.50) 



(3.51) 



which does not change the action. This transformation is more efficient compared 
to the heat-bath in terms of mobility, however, the move is not ergodic. Therefore 
we have to mix the over-relaxation update with the heat-bath update. 





(2,2) 



Figure 3.9: The {p, q) transformations in two dimensions 
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In the following sections we describe in detail the geometric transformations and 
the weights derived for the 5 {p, q) moves in four dimensions. For the transforma- 
tions where new field variables are created we derive also the heat-bath probability 
for the fields. To clarify the notation for the geometric transformations, consider 



first the 3 (p, g)-transformations in two dimensions depicted in figure 3.£. We use 
the notation tpip2P3 to denote a triangle which contains the vertices with the labels 
Pi,P2,P3- The first transformation 3) inserts a vertex into a triangle. We write 
this as 

t l23 ^ tl24 + tl34 + t234. (3.52) 

The inverse transformation Af(3^i) removes a vertex. With this procedure the num- 
ber of triangles grows (decreases) by 2, the number of links grows (decreases) by 3 
in the triangulation. 

The second transformation M^2,2) ^ips a link: 



tl34 + t234 



tl23 + tl24 



(3.53) 



This transformation is self-dual. Using the transformation at the 'same' location 
twice does not change the geometry. 

In the four dimensional case presented in the following sections we denote a four 
simplex spip2P3P4P5 by enumerating the labels of the vertices pi. 

3.5.1 Move 15, 51 

The move 5) inserts a vertex p6 into a four-simplex sl2345. This step generates 
4 four-simplices, 10 triangles, 5 links and, together with the links, five new field 
variables A{li6) = Xi,i = 1, . . . , 5. 



S12345 



S12346 + S12356 + sl2456 + sl3456 + s23456 
2 3 2 3 



(3.54) 





The change in the geometric part of the action is 

ASg = 4k4 - 10k2. 



(3.55) 



The order o{t) of the triangles t G sl2345 changes by The "old" fields therefore 
contribute 

(3.56) 



ASs{A) = ^{Aij + Ajk + Aki) 



The sum runs over all triangles which change the order by -|-1 in the transformation 

(3.57) 



-^(1,5)- The new fields contribute 



Sa[x] = 3 ^ {Aij + Xj - Xif 
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to the action. The sum runs over pairs i < j < 5 of the points of sl2345. For the 
probabihty H^p ^) (3.47) we have to calculate the volume V{x) ( |3.46 ) for the new 
fields. We find 



V 



-Sa[x] 



zero mode f 25 

With ( 3.47| ) we get for the Metropolis weight 



rk'^^'p U Yl (Aj + ^jk + Ak^f 

<ijk> 



n(i,5) = min <; 1, ^ !Jie exp [ J] {A, + A,k + A 

5 \ <ijk> 




(3.58) 



(3.59) 



for the geometrical part of the move. 

The inverse transformation removes a suitable vertex from the simplicial man- 
ifold. No new fields are created and hence Vf, = 1. The weight for this move 
is : 



11(5 1) = ( 1 



125 715 

— ——e ''exp 



(3.60) 



<ijk> 



After the geometric part of the transformation M^i^^-^ has been accepted the 
five new link-variables Xj = A{liQ) have to be assigned with a heat bath probability. 
One of them, say X5, may be set zero due to the gauge degree of freedom. The 
others are chosen with a Gaussian distribution : 



(3.61) 



P{xi,X2,X3,X4^,X5 = 0) ~ exp y-SY^i ji^jiAij + Xj - X. 

~ exp (— 3(2;-^Mx — 2w'^x)^ 



where 



M 



4 

-1 
-1 
-1 



-1 -1 

-1 -1 

4 -1 

-1 4 



(3.62) 



and 



Wi 



J2 ^^r (3-63) 

It is convenient to diagonalize the symmetric matrix M to decouple the fields. 
Denote the components of the vector in the eigenvector bases by yi,y2,y3,y4- The 
distribution for the in this bases is given by 



P{yi) ~ exp(-3(2/i - r2i)2) = i('u;i + i(;2 + i^s + ^^4) 

P(y2) ~ exp(-15(y2 - ^^2/5)^) = lim + W2 - W3 - W4) 

P{y^) ~ exp(-15(y3 - ^^3/5)2) = ^{wi - W2) 

P{yA) ~ exp(-15(y4 - ^^4/5)^) ^4: = ^{w^ - wa) 



(3.64) 
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This assigns a random value to the Ui with the distribution given above. The 
new fields Xi are obtained by a base-transformation : 

xi = ^{yi + y2 + V^ys) 
X2 = 
X3 = 
X4 = 



(yi + y2- V^ys) 
y2 + V^yi) 
y2 - V^Vi) 



1 

|(yi 



(3.65) 



3.5.2 Move 24, 42 

The move M(2,4) inserts one new link into the manifold. This generates also four 
triangles and two four-simplices: 

S1 3456 + S2 3456 ^ sl2345 + sl2346 + sl2356 + sl2456 (3.66) 





The change in the geometric part of the action is : 

ASg = 2k4 — Akappa2- (3.67) 

In the calculation of the volume V of the configuration space for the new field 
X = A{112) we use the notation ( p. 66 ): the new link 112 connects the vertices pi 
and p2, the simplices 513456 and s23456 have a common tetrahedron r3456. The 
order of the four triangles t G T3456 decrease by 1 in the move M(2,4) while the 
orders of the remaining 12 triangles in sl3456 and s23456 increase by 1. 

Denote by 

Pink = {Aii + + Akif (3.68) 



ijk — ~r -^jk "r -^ki) 

the value of the plaquette Pijt- One gets the following change in the action for the 
"old" fields : 



■i<i<j<k<6 2<i<j<6 

The volume of the state space for the new field is 



^ = yY^exp-3^(^i,-+^,-2 



i=3 

Collecting these results together one finds for the transition probability 



n 



(2,4) 



and 



n 



(4,2) 



min <! 1, ^ X y X g-^^f"'^^- 



min<!l,^xF-ixe^^«+^^^ 



(3.69) 

(3.70) 

(3.71) 
(3.72) 



for the inverse move M(4^2)- 

The heat-bath weight for the one new field is given by eq. ( 3.4^ ). 
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3.5.3 Move 33 



This move flips a triangle on the manifold. The transformation is self-dual and 
does not change the number of (sub)-simplices. Therefore we do not have to deal 
with a new field: 

V = 1. (3.73) 



S12456 + S13456 + s23456 ^ sl2345 + sl2346 + sl2356 (3.74) 



3 3 




In the transition probability one has to take into account the change A5s coming 
from the change of o{t). From |3.74| one can read off 

3 6 3 6 

i,j=l,i^j fc=4 j,k=l,jy^k i=4 

The action of the gravity sector is unchanged. Therefore the transition probability 
for the geometric part of the move is 

n(3,3) = min 1 1, ^ X e-^^» | (3.76) 
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Chapter 4 



Summary 



In this thesis we have analyzed two different models of random geometries, namely 
hyper-cubic random surfaces and four dimensional simplicial quantum gravity. Al- 
though these models are different, e.g. the local dimension of the models considered 
is two and four respectively, they have common properties. For example both mod- 
els possess a branched polymer phase. The appearance of this phase, however, is 
often referred to as a 'collapse to triviality' because the basically one dimensional 
structure does not describe extended d-dimensional {d = 2, 4) geometries in the 
continuum limit. One problem we addressed for both models was therefore to find 
modifications of the models with a different critical behavior. In both cases we find 
such modifications. 

The hyper-cubic random surface model is a discretization of bosonic strings 
embedded in a D-dimensional space. In three dimensions the model can also be 
interpreted as a fluid membrane. We have analyzed the phase diagram of the model 
with an extrinsic curvature term in the action. Numerically we find two phases 
in three and four dimensions: a weak coupling phase with a branched polymer 
behavior and a flat strong coupling phase. The Hausdorff dimension is d/^ = 4 in 
the branched polymer phase and dh = 2 in the flat phase. The transition is of 
first order, the critical coupling is finite. We have shown that the existence of this 
transition is not in contradiction with an earlier renormalization group analysis 
of the model, which under some assumptions indicated that the model is always 
in the branched polymer phase for finite coupling to the external curvature. One 
cannot take the continuum limit at the transition point because the transition is 
of first order. One can, however, try to smoothen this transition by introducing 
self-avoidance into the model. This is motivated by the observation that for the 
related model of Ising spins with a gonihedric action without self-avoidance of the 
dual surfaces evidence for a first order and with weak self-avoidance evidence for a 
second order phase transition was found |44, 45 1. 

We have solved an old question concerning the hyper-cubic random surface 
model. Earlier investigations indicated that the model with a certain local con- 
straint has a different critical behavior compared to the unrestricted model. We 
have shown that this is not the case. The true behavior of the constrained model 
is masked by very strong finite size effects. In the thermodynamical limit the con- 
strained and the unconstrained model have the same critical behavior. This shows 
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the universality of the model in the sense that the behavior of the model does not 
depend on details of the discretization. 

For the model of four dimensional simplicial quantum gravity we found that 
coupling more than two gauge vector fields to gravity suppresses the branched 
polymer phase. Our numerical data and a strong coupling expansion indicate that 
the branched polymer phase is replaced by a new phase of gravity. The Hausdorff 
dimension in this phase is close to four, the dimension of flat space. We find 
a negative entropy exponent 7 in this phase. This may indicate, similar to the 
model in two dimensions, that the whole phase is critical so that one can take 
a non-trivial continuum limit. One can hope to find gravitons in this model at 
the phase transition if the transition is of second order. However, the statistics of 
our numerical data does, at the moment, not allow to determine its order. Our 
ongoing numerical investigations of the model therefore concentrate on the nature 
of the phase transition. In the future the geometrical properties of the dominating 
manifolds at the transition point and in the new phase have to be analyzed. This 
will allow to answer the question of whether the dominating geometries describe, 
on large scales, a flat four-dimensional Euclidean space-time. If this is the case, we 
may have succeeded in finding a stable ground state of four dimensional euclidean 
quantum gravity. The observation that matter is required to stabilize the ground 
state of quantum gravity would then suggest that matter may be required for the 
very existence of a stable universe. 

A different question which we addressed for four dimensional quantum gravity 
concerns the role of topological excitations. It is not a priori clear whether topo- 
logical excitations should be taken into account in this model. We have checked 
numerically for three different topologies whether topological excitations are expo- 
nentially suppressed in the thcrmodynamical limit. We have shown that the leading 
contribution to the free energy does not depend on the topology and hence the topo- 
logical excitations survive in the thcrmodynamical limit. This result supports the 
conjecture that it is possible to define the model with coherent contributions from 
all topological excitations at a double scaling limit similar to the two dimensional 
model. 
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